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Abstract 

We study the Hamiltonian approach to 1 + 1 dimensional Yang-Mills theory in 
Coulomb gauge, considering both the pure Coulomb gauge and the gauge where 
in addition the remaining constant gauge field is restricted to the Cartan algebra. 
We evaluate the corresponding Faddeev-Popov determinants, resolve Gauss' law 
and derive the Hamiltonians, which differ in both gauges due to additional zero 
modes of the Faddeev-Popov kernel in the pure Coulomb gauge. By Gauss' law the 
zero modes of the Faddeev-Popov kernel constrain the physical wave functionals to 
zero colour charge states. We solve the Schrodinger equation in the pure Coulomb 
gauge and determine the vacuum wave functional. The gluon and ghost propagators 
and the static colour Coulomb potential are calculated in the first Gribov region 
as well as in the fundamental modular region, and Gribov copy effects are studied. 
We explicitly demonstrate that the Dyson-Schwinger equations do not specify the 
Gribov region while the propagators and vertices do depend on the Gribov region 
chosen. In this sense, the Dyson-Schwinger equations alone do not provide the full 
non-abelian quantum gauge theory, but subsidiary conditions must be required. 
Implications of Gribov copy effects for lattice calculations of the infrared behaviour 
of gauge-fixed propagators are discussed. We compute the ghost-gluon vertex and 
provide a sensible truncation of Dyson-Schwinger equations. Approximations of the 
variational approach to the 3+1 dimensional theory are checked by comparison to 
the 1 + 1 dimensional case. 
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1 Introduction 



In recent years there has been a renewed interest in Yang-Mills theory in the Coulomb 
gauge, both in the continuum approach [Tll2ll3]l4f5|6f7f8f9p0|n and on the 

lattice |17lll8|19|20j21j22|23| . Being non-covariant, Coulomb gauge is more cumbersome for 
perturbative calculations and, in fact, it has not yet been proved that Yang-Mills theory in 
the Coulomb gauge is perturbatively renormalisable. However, Coulomb gauge is superior 
to covariant gauges such as Landau gauge when it comes to nonperturbative investigations 
of the infrared sector of the theory. The reason is that in Coulomb gauge. Gauss' law can 
be explicitly resolved (and thus the separation of gauge dependent and gauge invariant 
degrees of freedom accomplished) [1]. This is a particular advantage in the Hamiltonian 
approach where the resolution of Gauss' law leads to a confining static potential between 
colour charges. A perturbative calculation of this potential allows the extraction of the 
running coupling constant [24|25| . Furthermore, in the so-called first order formalism of 
the functional integral approach, the resolution of Gauss' law in the Coulomb gauge leads 
to a cancellation of the Faddeev-Popov determinant and thus avoids in principle the in- 
troduction of ghost fields [T3|T4] . The Yang-Mills Hamiltonian in the Coulomb gauge is 
the starting point of a variational solution of the Yang-Mills Schrodinger equation in Refs. 
[6f7f8fT0] . Using Gaussian types of ansatze for the wave functional, minimisation of the 
vacuum energy density gives rise to a set of equations (similar to the Dyson-Schwinger 
equations) for various propagators and vertices. Imposing certain approximations, these 
equations have been solved analytically in the infrared |26j and numerically in the whole 
momentum regime [8|27fl2) . If the curvature of the space of gauge orbits expressed by the 
Faddeev-Popov determinant is properly included one indeed finds a linearly rising static 
potential and an infrared diverging gluon energy, both phenomena signalising confinement. 

In the derivation of the Dyson-Schwinger equations from the variational principle P|8] . 
a couple of assumptions and approximations are involved which are difficult to control 
in D = 3 + 1 dimensions. For this reason, we apply in the present paper the variational 
approach of Ref. [8] to Yang-Mills theory in D = 1 + 1, which can be solved exactly in the 
Coulomb gauge. With the exact solution at hand we then test various assumptions and 
approximations involved in the D = 3 + 1 dimensional case. 

Yang-Mills theory in 1 + 1 dimensions is trivial in fiat Minkowski space but becomes non- 
trivial on a compact manifold. We shall consider Yang-Mills theory on the space-time 
manifold 5*^x11 which is convenient for the Hamiltonian approach. On this manifold 
the Yang-Mills Schrodinger equation has been solved exactly with a (almost) complete 
gauge fixing f28]. This gauge fixing consists of the pure Coulomb gauge, which leaves in 
D = 1 + 1 a constant spatial gauge field, and an additional gauge condition, which exploits 
the residual global gauge invariance left in the pure Coulomb gauge to diagonalise the 
algebra-valued gauge field. We will refer to this (almost complete) gauge as the diagonal 
Coulomb gauge. The additional global gauge fixing is not implemented in the variational 
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approach in D = 3 + 1 dimensions. In D = 1 + 1, the pure and diagonal Coulomb gauges 
are here directly compared by identifying their Gribov regions and (previously unknown) 
boundary conditions on the wave functionals. The latter boundary conditions arise by 
paying proper attention to those zero modes of the respective Faddeev-Popov kernels that 
are not related to the field configurations on the Gribov horizon but to the incomplete 
gauge fixing. Thus, the exact vacuum wave functionals are derived in both gauges, and the 
propagators and vertices are calculated within the corresponding Gribov regions. These 
results for the Green functions are shown to be a particular solution of the exact Dyson- 
Schwinger equations defined on the first Gribov region. Other solutions to the Dyson- 
Schwinger equations are shown to exist and to be defined on a union of several Gribov 
regions, including many Gribov copies. We investigate which kind of approximation of the 
Dyson-Schwinger equations leads to which one of these solutions, for the sake of comparison 
to D = 3 + 1 where a truncation of the Dyson-Schwinger equations in mandatory. 

The calculation of the static colour Coulomb potential within a variational approach to 
D = 3 + 1 yields linear quark confinement [27] but includes certain approximations. We 
here calculate the static colour Coulomb potential in both the pure Coulomb gauge and 
the diagonal Coulomb gauge and compare the results for the Coulomb string tension to 
the gauge invariant string tension extracted from the Wilson loop. Along these lines, the 
quality of approximations in the D = 3 + 1 calculations can be estimated by comparison 
to the D = 1 + 1 case. 

The organisation of this paper is as follows: In section [2l we briefly review the Hamiltonian 
formulation of Yang-Mills theory in 1 + 1 dimensions on the space-time manifold x 
R, extract the physical conflguration space from the Wilson loop and identify the gauge 
invariant and gauge degrees of freedom. In section [3l we discuss incomplete gauge fixing 
in the Hamiltonian approach by the Faddeev-Popov method and calculate explicitly the 
Faddeev-Popov determinant of the pure Coulomb gauge, which differs from the one in 
the diagonal Coulomb gauge. In section we discuss the gauge fixing on in detail and 
determine the Gribov regions and the so-called fundamental modular region. The restriction 
of the configuration space to the fundamental modular region imposes boundary conditions 
on the wave functionals, which are also extracted in this section. In section [5l Gauss' law is 
explicitly resolved for both the pure and diagonal Coulomb gauges. The exact solution of the 
Yang-Mills Schrodinger equation in the pure Coulomb gauge is given in section [6l With the 
exact vacuum wave functional at hand, the exact propagators and vertices are calculated 
in section [71 In section [HI the potential between static colour charges is determined and 
the Coulomb string tension is extracted. In section [9], the Dyson-Schwinger equations are 
derived. Gribov copy effects on the Green functions are discussed in section [TOl The quality 
of a common truncation of Dyson-Schwinger equations is assessed in section [TH In section 
[T2l we apply the variational approach of Ref. [8] to Yang-Mills theory in D = 1 + 1 on 
5^ X R using the same variational wave functional in order to check the approximations 
made in Z) = 3 + 1. A short summary and some concluding remarks are given in section 
[T3l Some mathematical derivations are presented in appendices. 



3 



2 Hamiltonian formulation of 1 + 1 dimensional Yang— Mills theory on 5*^ x R 



We consider SU{Nc) Yang-Mills theory in 1 + 1 dimensions. In D = 1 + 1, the gauge field 
Afj_{x) is dimensionless while the gauge coupling constant g has dimension of inverse length. 
It is, however, convenient to absorb the gauge coupling g into the gauge field gA^ — » A^, 
so that A^ = A^j^T"- has dimension of inverse length. We will use antihermitian generators 
Ta, a = 1, . . . , Nl — 1 of the gauge group, satisfying 

\Ta-, Tb\ = fabcTc , (2.1) 

where fabc is the structure constant. For SU{2), the generators are related to the Pauli 
matrices Ta=i,2,3 by Ta = -^Ta. 

In 1+1 dimensional (fiat) Minkowski space Yang-Mills theory is trivial but becomes non- 
trivial on a compact manifold. The only compact space-time manifold with a canonical 
time is the cylinder 

(space) X R (time) . (2.2) 

On X R, the gauge invariant degrees of freedom are the spatial Wilson loops winding 
around S^, 

W[A] = ^trP exp |^-^dx^Ai(x°,x^)j , (2.3) 

which represent closed electric fiux lines. The spatial can be realized by consideringa 
finite interval on the x^-axis of length L and imposing the periodic boundary conditiorLU 

A^{x^,x' = L)=A^{x'',x' = 0). (2.4) 

In the absence of fermions in the fundamental representation this boundary condition 
remains intact under gauge transformations 

A^^ = UA^U^ + Ud^U^ (2.5) 

satisfying the boundary condition 

f/(x°, x^ = L) = Zn{x^)U{x", = 0) , (2.6) 

where Zn, n = 0, 1, . . . , Ai'c — 1 is an element of the centre Z{Nc) of the gauge group. 



Throughout the paper we will work in the canonical Hamiltonian approach [30] and impose 
the Weyl gauge 

^0 = 0. (2.7) 



^ The periodic boundary condition can be taken without loss of generality since the SU{Nc) 
bundle over 5^ is trivial as noticed in Ref. 1291. 
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The gauge transformation required to bring a periodic gauge field y4o(x°,x^ = L) 
Aq{x^,x^ = 0) into the Weyl gauge, 



U\x'^,x^) = P 



exp 



/ 



V 



dtAo(t,x^] 



(2i 



is also periodic and thus within the class l\2M . On the flat space R, an analogous gauge 
transformation 



V\x'^,x^) =Pexp - / dsAi{x\s) 



(2.9) 



could also gauge away the field Ai, i.e. AY = 0. However, the gauge transformation V (12. 9p 
is not within the class (12.61) and therefore not allowed on S^. It is the compactification of 
space R — 5> S"^ which makes the theory non-trivial. 



The canonical quantisation is carried out at a fixed time x", so that Ai{x'^ — Hxed, x^) is the 
only field "coordinate" left. To simplify the notation, we will omit henceforth the spatial 
index i = 1 and write x^ — x , di ^ d , A"(x^) — * A"-{x) etc. The Hamiltonian of 1 + 1 
dimensional Yang-Mills theory then reads 

2 

J dxW{x)W{x) , (2.10) 

where 

U'^^x) = ---^ (2.11) 

is the momentum operator, which represents the electric field. Note in 1+1 dimensions, 
there is no magnetic field and hence no potential term in the Hamiltonian. 

Having quantised the theory in the Weyl gauge, the classical Gauss law DU = pis lost from 
the quantum equations of motion. Enforcing Gauss' law as an operator identity contradicts 
the canonical commutation relations!!] One therefore imposes Gauss' law as a constraint 
on the wave functional "^{A) 

D{x)U{x)^{A) = p{x)'^{A) . (2.12) 

Here, 

D{x) = d + A{x), A{x) = A''{x)fa, (2.13) 



^ Alternatively, the Dirac bracket formalism can be used to quantise the theory after fixing the 
gauge transformations generated by Gauss' law. The Gauss law can then be imposed as an operator 
identity. This was shown in 3+1 dimensional Minkowski space to lead to the same energy spectrum 

mi. 
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denotes the covariant derivative in the adjoint representation with 

{fa)bc = fbac (2.14) 

being the generators in the adjoint representation. Furthermore, p{x) denotes the colour 
density of matter fields (or external sources). The operator on the l.h.s. of Eq. (I2.12p . DU, is 
the generator of so-called "small" gauge transformations (see section HI) and in the absence 
of matter fields, p(x) = 0, Gauss' law forces the wave functionals ^^(^4) to be invariant 
under the small gauge transformations U, 

^{A^) = ^{A) . (2.15) 



Instead of working with gauge invariant wave functionals, it is more convenient to explicitly 
resolve Gauss law by fixing the gauge [1] and for this purpose the Coulomb gauge 

dA{x) = (2.16) 

is particularly convenient. In one spatial dimension, the Coulomb gauge constrains the 
gauge field to spatially constant modes and the field theory reduces to quantum mechanics 
in these constant modes. The gauge transformation Q which brings a given gauge field A 
into the Coulomb gauge dA^ = can be chosen 

n^x) = V\x){V{L))t: , (2.17) 

where V{x) is defined by Eq. (12. 9p with the x^-dependence suppressed. Contrary to V{x), 
Q{x) is periodic for periodic A{x), and thus within the class of allowed gauge transforma- 
tions (12.611 . The gauge-transformed field 

A^ = ylnV{L) (2.18) 

is space-independent and thus obviously satisfies the Coulomb gauge condition. 

The Coulomb gauge (12.16^ (together with Weyl gauge (I2.7p ) still leaves invariance with 
respect to global (space and time independent) gauge transformations, which can be ex- 
ploited to diagonalise the constant gauge field 

A = A'^Ta ^ A'^^Ta, = A^^ag ■ (2.19) 

Here T^y, = 1, N^ — l denotes the generators of the Cartan sub-algebra which 
can be chosen to be diagonal. Equation ( 12.1911 fixes the global transformation [/ up to a 
constant element of the Cartan sub-group, i.e. even after implementing the gauge condition 
A = diagonal in addition to the pure Coulomb gauge, there still is a residual global abelian 
f/(l)^'=~^ gauge symmetry. Moreover, there still is a discrete symmetry left: The above 
gauge conditions (12.16P and (12.191) do not fix the so-called Weyl symmetry consisting of 
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the NJ permutations of the Nc eigenvalues of A°'Ta. In the case of SU{2)^ where the two 
eigenvalues have the same modulus but opposite signs, the Weyl symmetry switches the 
signs of the eigenvalues. 

Note that we can re-express the two gauge conditions f l2.16p and (I2.19p as 

= = 0, (2.20a) 

/""[A] = M"«(x) =0, (2.20b) 

i.e. the Coulomb gauge condition is imposed only on the abelian part A"-° while the non- 
abelian components of the gauge field denoted by the index a ^ Qq vanish. In the following 
we will refer to this gauge as diagonal Coulomb gauge, while Eq. (12.161) will be called pure 
Coulomb gauge. 

To simplify the explicit calculations, we will confine ourselves henceforth to the gauge 
group SU (2) where the structure constants coincide with the totally anti-symmetric tensor 
jafec _ ^abc^ ^ _ 2^ 3^ and we will choose the generator of the Cartan sub-group to be 
given by oq = 3. For SU{2), the A"" can be interpreted as the Cartesian components of a 
vector A = A^e^ in a 3-dimensional Euclidean space with Cartesian unit vectors e^^i 2,3. 
In the diagonal Coulomb gauge specified by Eq. (I2.20p the colour vector A is rotated into 
the direction of the positive or negative 3-axis 

A = A^es , A^ = ±|A| . (2.21) 

The two signs differ by a Weyl refiection A —A. Thus the gauge transformation from 
the pure Coulomb gauge into the diagonal Coulomb gauge is given by rotating the colour 
vector A into the positive 3-direction 63, possibly followed by a Weyl refiection A^ —A^. 

The modulus |A| of the constant Coulomb gauge field represents the only gauge invariant 
degree of freedom. In fact, the only physical observable of the theory, the spatial Wilson 
loop W winding (non-trivially) around the whole space manifold S^, Eq. (12. 3p . is easily 
calculated in (both pure and diagonal) SU{2) Coulomb gauge to be given by 

W =]-iY exp(-AL) = cos?9 , (2.22) 

where we have introduced the dimensionless variable 

= ^ , (2.23) 

for later convenience. W attains all possible values in [—1, 1] when d traverses the interval 

< ^9 < vr . (2.24) 

Equation (12.24^ defines the physical configuration space of the theory. We will later recover 
the interval (I2.24p as the so-called fundamental modular region. 
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To separate the constant Coulomb gauge field A into (global) gauge invariant andgauge 
dependent parts, it is convenient to use the spherical coordinates |A|,6', to writqiJ 

A=\A\k{d,(P) = Y^k{e,(P), (2.25) 

where 

A(6', (j)) = sin 9 (cos (pei + sin (p 62) + cos 9 63 (2.26) 
is the radial unit vector. 

The primary aim of the present paper is to use 1+1 dimensional Yang-Mills theory as 
testing ground for the variational approach developed in 3 + 1 dimensions [8], assessing the 
approximations introduced there. In studies of Coulomb gauge Yang-Mills theory in higher 
dimensions, merely the pure Coulomb gauge is fixed. Therefore, we will here mainly focus 
on the pure Coulomb gauge as well. Moreover, the diagonal Coulomb gauge will be used to 
investigate the effect of a complete gauge fixing on the Green functions of the theory. 



3 Gauge fixing in the Hamiltonian approach by the Faddeev-Popov method 



In the Hamiltonian approach in Weyl gauge, there is no need to fix the residual time- 
independent gauge symmetry. In principal, one can work (in the absence of external colour 
charges) with gauge invariant wave functionals, which trivially satisfy Gauss law. It is only 
a matter of convenience that one prefers to fix the gauge. Furthermore, as will be explicitly 
shown in the context of the resolution of Gauss' law (see section [5]), the gauge fixing needs 
not to be complete, i.e. any partial gauge fixing is allowed in the Hamiltonian approach. In 
any case, the wave functionals have to be invariant under the residual gauge symmetries 
unfixed by the gauge condition. In the pure Coulomb gauge, the wave functionals have to 
be invariant under global colour rotations U = const, which are not fixed by that gauge. 

Gauge fixing is accomplished in the Hamiltonian approach by applying the Faddeev-Popov 
method to the functional integral over the spatial gauge fields defining the scalar product 
in the Hilbert space of wave functionals. Consider the matrix element of a gauge invariant 
observable 0[A^] = 0[A] 

(^i|C>[A]|^2) = JVA<I/1[A]0[A]<I/2[A]. (3.1) 
The Faddeev-Popov method amounts to inserting into the functional integral the identity 

1 = Bet M[A] J Dfi{Q) 6{r[A^]) (3.2) 

^ For a 3- vector A, the caret denotes as usual the unit vector A = A/|A|, while for algebra- 
and group-valued quantities the caret means the adjoint representation. 
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where ri = e , 9 = G'^T^, denotes a gauge transformation and 

-'(-'^^ 

is the Faddeev-Popov kernel. Furthermore, for a complete gauge fixing, Dfi{Q) denotes the 
Haar measure of the gauge group. Inserting the identity (13.21) into Eq. (13. ip and exploiting 
the gauge invariance of both the wave functional and the observable, one finds after a 
change of the integration variable 

(*i|0[A]|^2) = / VADetM[A] 5{r[A])^i*[A] 0[A] ^2[A] J Dfi{Q) (3.4) 

where the integration of the gauge group is now explicitly separated, yielding a (infinite) 
constant which can be absorbed into the normalisation of the wave functional. 

If the gauge condition /"[A*^] = does not fix the gauge completely, there are directions 
in the space of gauge transformations, along which the gauge-fixing functional /"^[A^] does 
not change, i.e. Sf"'[A^]/6Q^ = 0, and the tangent vectors corresponding to these directions 
represent zero modes of the Faddeev-Popov kernel (13. 3p . In order that the identity (13. 2p 
holds and thus the Faddeev-Popov method works one has to exclude these zero modes and 
integrate only over the subspace of gauge transformations which are fixed by the gauge 
condition /"[A^] = 0. Furthermore, the gauge fixing is defined only in the region of gauge 
field configurations A where the Faddeev-Popov determinant Det A^[y4] is non-zero. This 
will be important for the application of the Faddeev-Popov method given below. 

The explicit calculation of the Faddeev-Popov kernel is most easily accomplished by notic- 
ing that the Gauss law operator DH is the generator of (so-called smalQ) gauge transfor- 
mations 

r[A'']=g{Q)r[A]g-\Q), (3.5) 
^(6) = exp (^-i J dx (]D''\x)Q\x)) W{x) 

For infinitesimal gauge transformations Q{5Q) = 1 + dx 5Q°'{x){DIlY{x), we have ac- 
cordingly 

r[^^](a;) = nA]{x) + j dySQ\y) [t{DIl)\y), r[A]ix)] (3.6) 

and comparison with the Taylor expansion of /"[A^] in powers of 6Q°'{x) reveals that the 
Faddeev-Popov kernel ( 13. 3p can be expressed as 

M''\x,y) = ib'\y) [ff(y),r[A](x)l . (3.7) 



^ See section m 
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Below, we use this relation to determine the Faddeev-Popov kernels for the two gauges 
considered above and calculate the corresponding Faddeev-Popov determinants. 



3.1 Diagonal Coulomb gauge 



Applying the relation (13 .7^ to the diagonal Coulomb gauge condition f l2.20p . we obtain the 
Faddeev-Popov kernel (6 = 1,2) 



-b'''\x)d''5{x,y) 



(3.8a) 
(3.8b) 



On the gauge shell (where the gauge condition is fulfilled), the Faddeev-Popov kernel 
becomes block diagonal and reads in the Cartesian basis 



M{x,y) 



( 



V 



-dx 









-dx 














\ 



S{x,y) 



i 
idx 



iD{x) 5{x, y) 



(3.9) 



The Faddeev-Popov determinant becomes 



T)eiM = Det(i) {-idy^ Det {iD 



(3.10) 



The first factor is a (divergent) irrelevant constant, which drops out from the expectation 
values and can be absorbed into the functional integral measure, see below. Note, this 
determinant is defined in the (one-dimensional) abelian colour sub-space a = = 3 only, 
indicated by the sub-script (1). To evaluate the second factor Det (iD), we consider the 
eigenvalue equation 

iD''''[A]{x)ifi\x) = \ifi''{x) (3.11) 

From the representation (13.31) . it is clear that the eigenfunctions of the Faddeev-Popov 
kernel have to satisfy the same boundary condition as the gauge angles 0"(x). Since the 
gauge transformations on the spatial have to satisfy periodic boundary conditions 6(x + 
L) = Q{x), the eigenfunctions of the Faddeev-Popov kernel have to be periodic as well. 
Since d and D commute on the gauge-fixed manifold, the operators d, D and have 
common eigenfunctions and we can impose the periodic boundary condition also on the 
eigenfunctions of zD, 

(/?"(x + L) = ¥;"(x) . (3.12) 

With this boundary condition, which does not mix the different colour components, and 
with the fact that the gauge field is spatially constant, the eigenfunctions of iD factorise 
in space and colour dependent parts. Using the results of appendix [Al and Eq. (I2.2ip . the 
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eigenfunctions are given by 



(fn,ai^) = {x,a\n,a) = {x\n){a\a) (3.13) 

where {a\a) = e^ are (the Cartesian components of) the polar unit vectors e^r (IA.15I1 and 

1 Oirv 
{x\n) = -= e"'^"" , fc, = -— , n e Z (3.14) 
V L L 

are plane waves (periodic in L). The corresponding eigenvalues read 

An,a = kn + (tA^ . (3.15) 

The two eigenvalues A„,cr=±i correspond to the non-abelian block (a = 1, 2) in the upper left 
corner in the Faddeev-Popov matrix (13.91) while the eigenvalue Xn,a=o = corresponds to 
the abelian colour direction oq = 3 (see Eq. (IA.15I1 ). Furthermore, the zero mode ipn=o,a=o = 
eo-=o/V^ = e^/VL (A„=o,o-=o = 0) represents the tangent vector (to the gauge orbit) 
corresponding to the infinitesimal global U{1) colour rotation (O^ = const) which is not 
fixed by the gauge condition (12.201) . As discussed above, this mode has to be excluded from 
the spectrum of the Faddeev-Popov kernel (13. 9p . whose eigenvalues are given by 

, iKi a , cr = ±1 

K.= { _ / _ • (3.16) 

Note also that the eigenmodes (pn=o,a=±i corresponding to the global SU{2)/U{1) gauge 
transformations which are fixed by the gauge condition (12.20p do not give rise to zero 
eigenvalues, A„=o,o-=±i = 7^ 0. 

Excluding the zero mode n = a = (indicated in the following by a prime), we obtain 
with = ±|A| (see Eq. ([221])) 

oo oo 

Det' {iD[A])= n ' n A„,,=(n^™) n 

n=— oo a=0,±l m^O n=— oo 



nfc^)iAr(n(^n-iAi 

rrif^O n=l 



2 



oo I A |2 ^ 2 



( n *y ( iAi n (1 - ^) ) (3.17) 



m^O \ n=l "'n 

The first factor represents Det (it?) with the zero mode n = excluded. Using 



sinx = X TT fl - f — ) ) (3.18) 

n=i ^ ^Tin^ ^ 
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we obtain 



Det' iiD[A]) = Det (td) [jf sin^ (^^^ = Det (td) (|-)' sin^ ^ . (3.19) 

All field-independent factors in the Faddeev-Popov determinant (13.101) can be absorbed in 
the functional integral measure. We thus arrive at the Faddeev-Popov determinant J'd of 
the diagonal Coulomb gauge (I2.20p 

Jd '■= 7 — -oNco — , , o = sm {} . (3.20) 

Det(i) (-2(9)33 Det(29) (|)2 ^ ' 

Let us stress that it was absolutely crucial to exclude the gauge modes which are not fixed 
by the gauge condition. Otherwise the Faddeev-Popov determinant would have vanished 
identically. 

The Faddeev-Popov determinant J'd has zeros at (note that by definition ^ > 0) 

■i? = ^^ = n7r, n = 0,1,2,... (3.21) 



and thus divides up the gauge-fixed configuration space into regions where the Faddeev- 
Popov method of gauge fixing is defined. These so-called Gribov regions are given by 

nvr < (n+ l)7r . (3.22) 

Recall that in the diagonal Coulomb gauge = ±|A|, so that in this variable the Gribov 
regions are given by 



A, 



2nn , , , 2(n + l)7rl , , 

< |^3|< L \ ^ n = 0,l,2,... . (3.23) 



The boundaries of the Gribov regions, the Gribov horizons, are given by the discrete mo- 
menta kn (I3.14p . We will return to the discussion of the Gribov regions in section [H 

Note that the Faddeev-Popov determinant J'd f l3.20p vanishes also at the (classical) per- 
turbative vacuum A = 0. This is not surprising since the diagonalisation of the gauge field, 
and thus the diagonal Coulomb gauge, is ill-defined for A = 0. Therefore, this gauge is not 
suitable for perturbation theory. 
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3.2 Pure Coulomb gauge 



Let us now consider the pure Coulomb gauge (12.161) which leaves the global gauge trans- 
formations unfixed. The Faddeev-Popov kernel (13.70 is then given by 

M''\x, y) = {-b'^'lA] {x)d^) 5{x, y) . (3.24) 

Since the gauge field in the pure Coulomb gauge is related to the one in the diagonal 
Coulomb gauge by a global gauge transformation, we can express the covariant derivative 
Z)'^''[A] in the pure Coulomb gauge by the one in the diagonal Coulomb gauge. Let 
denote the global gauge transformation which rotates the colour vector A into the positive 
3-direction, i.e. 

A'^fa = lApf^U^ . (3.25) 

Then we have (see appendix [All 

iD[A] = UiD[\A\T3] (3.26) 

where U is the adjoint representation of U. Here, i-D[|A|T3] is the covariant derivative in 
the diagonal Coulomb gauge (with A^ = ±|A|) whose eigenvalues and eigenfunctions were 
determined in the previous subsection. From Eq. ( I3.26p . it follows that iID[A°'Ta] has the 
same eigenvalues A„o- (I3.15P as zZ)[|A|T3], i.e. 

K,a = K + (r\A\ , (3.27) 

and that the eigenfunctions \ip) of zZ)[74] are related to the eigenfunctions \ip) (13.131) of 
iD[\A\T3] by 

^l^{x) = Uab^lA^) = {x\n)u: , (3.28) 
where {x\n) are the periodic plane waves (13.141) and we have defined (see appendix [All 

u-:={a\U\a) = {a\U\b){b\a) = Ua,e'^ 

= {a\T){T\U\a) = eMA<l^,0,0). (3.29) 

Since «i)[y4] ( I3.26p has the same eigenvalues as ?D[|A|T3], and since det U = 1, one would 
expect that the Faddeev-Popov determinant in the pure Coulomb gauge gauge is, up to 
an irrelevant constant factor Det (i(?)/Det(i) {—id)^^, the same as in the diagonal Coulomb 
gauge (|2.20p considered above. However, since the pure Coulomb gauge does not fix the 
global gauge transformation U, we have to exclude the constant eigenmodes from the 
Faddeev-Popov kernel. These are given by the eigenfunctions with n = and all a. In 
addition to the zero mode n = 0, a = excluded already from the Faddeev-Popov kernel 
of the diagonal Coulomb gauge, one has to exclude here also the gauge modes n = 0, a = 
±1, corresponding to the non-zero eigenvalues A„=o,ct=±i = =1=1^1. Although these modes 
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(a) (b) 

Fig. 1. The first Gribov region for (a) the diagonal Coulomb gauge and (b) the pure Coulomb 
gauge. 

correspond to non-zero eigenvalues of i-D[A] they give rise to zero modes of the full Faddeev- 
Popov kernel of the pure Coulomb gauge, 

{-D[A]d)(^n,a = K,a^n.p , ^n.p = K.aK (3.30) 

since A„=o,(7 = 0. Note that these zero modes form precisely the global gauge transformation 
U which diagonalises the constant Coulomb gauge field A = A°-Ta (cf. Eqs. ( 13.26p . (lA.lf l). 

^l=,A^) = ^(«|f>l^) • (3.31) 

Omitting these modes in calculating Det {—Dd), equivalently to Eq. (13.171) . we find for the 
Faddeev-Popov determinant jTp in the pure Coulomb gauge 



Det{-Dd) _ sin' (^) _ sinH 



Det(-a2) " (^l^y - ^2 • (3-32) 

This determinant differs from the one of the diagonal Coulomb gauge f l2.20p by the denom- 
inator (cf. Eq. (|3.20p ) and does not vanish at the perturbative vacuum A = 0. Furthermore, 
since \im\A\^o dj^p / d\A\ = the space of gauge orbits is flat near the perturbative vacuum 
in this gauge. The first zero of the Faddeev-Popov determinant occurs at -d = \A\L/2 = vr 
which defines the first Gribov horizon to coincide (in the variable with the one in the 
diagonal Coulomb gauge (I3.22p . However, in the pure Coulomb gauge the direction of the 
gauge field A is not fixed, so the Gribov regions Qn are given here by the spherical shells 



n„. := A 



[n-1)— <\A\<n — \ , n = 1,2,3,.... (3.33) 

In particular, the first Gribov region Vti is given by the sphere |A| < The first Gribov 
regions of both the diagonal and the pure Coulomb gauge are illustrated in Fig. [H The 
Gribov regions will be discussed in more detail in the upcoming section. 
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4 Gribov regions and boundary conditions on the wave functionals 

In this section, the configuration spaces of both the pure and diagonal Coulomb gauges are 
examined. (For valuable discussions on this topic, also see, e.g., Ref. [32] and references 
therein.) The topology of gauge transformations relating the various Gribov regions of the 
pure Coulomb gauge and the diagonal Coulomb gauge is discussed here, the fundamental 
modular region is specified, and boundary conditions on the wave functionals are given. 



4.1 Topology of gauge transformations on S*^ 

Since the first homotopy group Ili{SU{Nc)) is trivial, there are no "large" gauge trans- 
formations in one (compact) spatial dimension, i.e. on S^. However, in the absence of 
matter fields in the fundamental representation, the gauge group is in fact SU{Nc)/Z{Nc). 
This is because the gauge field, living in the adjoint representation, is invariant under 
centre gauge transformation^ Z{x) e Z{Nc). Since Ui{SU{Nc)/Z{Nc)) = Z{N^), there 
are topological non-trivial gauge transformations falling in Nc different topological classes 
and consequently there are N^. distinct classical vacuum configurations Af^n) = ^ri)'9f^(n) 
between which quantum tunnelling occurs. In the fundamental representation, the f/(„) sat- 
isfy the boundary condition (12. 6p , with x° fixed, and are specified by the Nc centre elements 
Zn G Z{Nc). The U(^n=o), belonging to the trivial centre element Zq = 1 are periodic and 
form the "small" gauge transformations. The remaining ones are "large" transformations. 
Since (Zn)'^" = 1, Nc successive large gauge transformations U(n) belonging to the same 
Zn yield a small gauge transformation. The fundamental representation of SU{Nc)/Z{Nc) 
is the adjoint representation of SU{Nc), which we will denote in the following by a caret. 
The adjoint representation U = exp(9aTa) is related to the fundamental representation 
U = exp(9aTa) (with the same 9^ !) by 

U^TaU = UabTt , Uab = "2 tr {u^TaUn) . (4.1) 

From this representation it is explicitly seen that the adjoint representation U is cen- 
tre blind. Accordingly, the allowed gauge transformations satisfying (in the fundamental 



^ Since the centre Z{Nc) is a discrete set, centre gauge transformations Z{x) have to be piece- 
wise constant. At the jumps of Z{x) (from one centre element to another) the inhomogeneous 
term of the gauge transformation, Z{x)dZ^x), represents so-called ideal centre vortices |33j34| . 
The centre vortices of the fundamental representation, Z{x)dZ^x) = C"'{x)Ta, become invis- 
ible Dirac sheets (strings) in the adjoint representation Z{x)dZ\x). In fact, a centre element 
Z{x) = exp(— /i'^(x)Ta) G Z{Nc) (in the fundamental representation, with //"(x) being a co-weight 
vector) becomes the unit matrix in the adjoint representation, Z{x) = exp(— /i°'(x)Ta) = 1, so 
that the inhomogeneous term Z{x)dZ^ (x) disappears in the adjoint representation. 
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representation of SU{Nc)) the boundary condition (12.61) are periodic in the adjoint repre- 
sentation 

U{x + L) = U{x) . (4.2) 

To be more specific consider the gauge group SU{2) whose centre is Z{2) = {1, —1}. The 
allowed gauge transformations satisfying Eq. (12. 6p are either periodic (belonging to Zq = 1) 
or anti-periodic (belonging to Zi = — 1). In the absence of matter fields in the fundamental 
representation, the true gauge group is SU{2)/Z{2) = 5*0(3). Since Ili{S0{3)) = Z{2)^ 
there are two inequivalent sets of gauge transformations allowed by Eq. ( 12.611 . The periodic 
ones U{Q){L) = f/o(0) belonging to the trivial centre element Zq = 1, form the small 
gauge transformations, which can be smoothly deformed to unity. The anti-periodic gauge 
transformations, corresponding to the non-trivial centre element Zi = —1, cannot be 
smoothly deformed to unity and are called large. To provide some explicit examples consider 
the gauge transformation 

U{x) = exp {uj{x)n ■ T) = exp (—i — h-Tj = cos — — i n • r sin — - (4.3) 

with some constant unit vector fi. A small (periodic) gauge transformation is obtained 
when ijj{x) satisfies the boundary condition 

a;(x + L) = u;(a;) + 4n7r , (4.4) 

while a large (anti-periodic) gauge transformation follows for 

u;(x + L) = cj(x) + (2n + l)27r . (4.5) 

The corresponding adjoint representations {iT,^hc = Cfeac) 

{U{x))ah = (e'^"''^)^^ = Sab COS uj + nanb{l -cosuj) + €^06^^ siu w (4.6) 
are periodic in both cases. 



4-2 Gribov and fundamental modular regions 



In section [3l we found that in the pure Coulomb gauge the Faddeev-Popov determinant 
has zeros at ?9 = -^^^ = nn. Hence, in the 3-dimensional space of constant gauge orbits 
the Gribov horizons dVtn are given by the surfaces of spheres around the origin with radius 
|A| = ^(^? = nvr), see Eq. (13:331 ). Note also that in the limit L ^ cx) the first Gribov 
region Vti shrinks to the point A = 0, in agreement with the fact that 1 + 1 dimensional 
Yang-Mills theory becomes trivial on a flat space-time manifold. 

The pure Coulomb gauge is not a complete gauge fixing since it leaves invariance with 
respect to global gauge transformations, which form the zero modes of the Faddeev-Popov 
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Fig. 2. Illustration of the displacements of a field configuration A from inside the first Gribov 
region to Gribov copies in the neighbouring Gribov regions by the large gauge transformations 

y±i = (y(A))±i. 

kernel. However, the global gauge transformations are not the only symmetries left. In the 
following we will carefully examine the residual symmetries left after pure Coulomb gauge 
fixing. Due to the existence of these residual symmetries, the first Gribov region Qi cannot 
yet be the fundamental modular region which, by definition, contains only a single copy 
of each gauge orbit. Since the pure Coulomb gauge is contained in the diagonal Coulomb 
gauge, the fundamental modular region is the same in both gauges, while the Gribov regions 
are, of course, different. The determination of the fundamental modular region is necessary 
in order to identify the symmetry relations (i.e. boundary conditions) to be fulfilled by the 
wave functionals. 

Consider the large (i.e. anti-periodic) gauge transformation, cf. Eqs. f l4.3p and (14. 5p . 

V{A) = exp (^^x A ■ ' ^ " ^ ^'^''^^ 

which shifts a constant gauge field A = |A|n along its direction n in colour space by 
multiples of 

A ^ A^ = A - ^A = (|A| - ^) A . (4.8) 

Obviously, the transformed configuration still satisfies the (pure or diagonal) Coulomb 
gauge condition if the original does. 



The large gauge transformation (14.81) maps a configuration from the n*'^ Gribov region 
to a configuration (on the same ray through the origin in colour space but) within fl„--i 1^1 
In this way, any gauge configuration in one Gribov region has a unique copy in every other 
Gribov region and all Gribov regions are homeomorphic to each other. Fig. [2] illustrates the 
shifting of a particular gauge configuration from the first Gribov region to the neighbouring 
ones. Furthermore, the large gauge transformation (14. 8p maps a configuration of the n*'^ 
Gribov horizon dVL^ to the configuration on the same ray through the origin on dVLn-i- 

^ For a peculiarity of n = 1, see below. 
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In particular, the configurations on the first Gribov horizon dVli, where |A| = are 
mapped to the vacuum A = 0. This shows that all configurations of (all) Gribov horizons 
are equivalent under large gauge transformation to the vacuum A = 0. In particular, since 
two successive large gauge transformations form a small one, all configurations on a given 
Gribov horizon are related by small gauge transformations and are thus gauge copies of 
each other. Let us stress, however, that the configurations on the first Gribov horizon 
dfli (and on all dfln with odd n) are not equivalent to the vacuum A = with respect 
to small gauge transformations, they are related to A = by a non-trivial (large) gauge 
transformation. In the quantum theory tunnelling between these vacua will occur. This 
tunnelling will be entirely accounted for by solving the Schrodinger equation. 

By a large gauge transformation (14.71) . a configuration A inside the first Gribov region with 
|A| < ^ is also mapped to a copy A^ with J < |A^| < A^ = —A and vice versa. One 
could eliminate the large gauge transformation by restricting the modulus of the gauge 
field to 

0<|A|<y. (4.9) 

However, Gauss' law only enforces the wave functionals to be invariant under small gauge 
transformations while they can transform according to an arbitrary representation of the 
symmetry group under large gauge transformations. For example, in the colour singlet 
sector under a large gauge transformation the wave functional can acquire a non-trivial 
phase e*", which is well-known from the 0-vacuum in D = 3 + 1. Therefore we will not 
remove the large gauge symmetryQ 

The first Gribov region in the pure Coulomb gauge, fii, is a ball around the origin 
with radius |A| = bounded by the first Gribov horizon dfli, which is the 5*2. Since all 
configurations on this Gribov horizon are equivalent with respect to small gauge transfor- 
mations, we have to identify all points of the first Gribov horizon, 5*2, which compactifies 
the first Gribov region to S3, which is the manifold of the SU{2) group. This shows that the 
configuration space of Yang-Mills theory on S"^ x IR is the gauge group manifold itself. The 
mapping from the configuration space into the gauge group is provided by the (untraced) 
spatial Wilson loop winding around S^. In the fundamental representation we have 

W[A] = Pexp 1^ j^dxA^ = exp(LA • T) = 1 cost? + 2A ■ T sin ?9 , (4.10) 

. |A|L , . , 27r 

<' = ^, 0<|A|<-. 



^ If one restricted the configuration space to |A| < ^, from the large gauge transformations only 
a residual discrete symmetry on the new Gribov horizon |A| = ^ would be left, which is given 
by the displacement transformations V{A) (|4.7p and which relates the antipodal points ±f A. 
Identifying these antipodal points, which are equivalent by the displacement transformation (14. Sh 
the Gribov region, the ball with radius ^, becomes the group manifold of 50(3), which is the 
gauge group in the absence of fundamental charges. 
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The Gribov horizons dfln with odd n (and in particular dfli) are mapped onto the non- 
trivial centre element 

ly |A| = (2n + 1)^1 = -1 . (4.11) 

In higher dimensions, field configurations A, non-trivially linked to a closed loop C for 
which the corresponding Wilson loop Vr[yl](C) equals a non-trivial centre element, are 
referred to as centre vortices. In this spirit, the field configurations on the Gribov horizons 
dfln with odd n represent centre vortices, in agreement with the general observation that 
centre vortices are on the Gribov horizon p3j, due to their larger symmetry. 



The first Gribov region Qi of the pure Coulomb gauge can be restricted further by im- 
plementing the diagonal Coulomb gauge (12.211) . In the diagonal Coulomb gauge, the first 
Gribov region is given by 

A^L 

-7r<^<7r, (4.12) 

see Eq. (13.231) . and is obviously a subset of fli. There are still gauge copies within the first 
Gribov region of the diagonal Coulomb gauge, due to the fact that there remains a residual 
discrete gauge symmetry, the Weyl reflection 

Removing this symmetry by identifying configurations of opposite sign, 

A^ = \A\ = ^, (4.14) 

reduces the configuration space (14.121) to the fundamental modular region 

27rl 



A-" e 



t9G[0,7r]. (4.15) 



This physical configuration space was already found from the gauge-invariant spatial Wil- 
son loop, see Eq. (12.241) . and reduces the configuration space to the genuinely gauge in- 
variant degree of freedom -d = |A|L/2. The Faddeev-Popov determinant JId{^), see Eq. 
( 13.20p . is gauge invariant as well, since it corresponds to the Faddeev-Popov determinant 
of a completely fixed gauge. Such a Faddeev-Popov determinant is gauge invariant due to 
the invariance of the Haar measure P/i(f/). 

The gauge condition that immediately rotates the constant colour vector A of the pure 
Coulomb gauge into the positive 3-direction is accomplished by the gauge transformation 
U given by Eq. (lA.QD . In this gauge the Faddeev-Popov determinant has zeros (Gribov 
horizons) at {} = nvr , n G N and the Gribov regions are given by the one-dimensional 
intervals 

mr<'d < {n + l)n , n = 0,l,2,... . (4.16) 
For such a gauge, the first Gribov region, n = 0, coincides with the fundamental modular 
region (14.151) . 
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4-3 Boundary condition on the wave Junctionals 



Let us now discuss the implications of the residual gauge symmetries on the wave func- 
tionals. In general, by Gauss' law the residual gauge symmetries which correspond to small 
gauge transformations that are not fixed by the gauge considered have to be respected by 
the wave functional. In the pure Coulomb gauge global gauge invariance is left unfixed and 
consequently the wave functional has to respect this symmetry, i.e. 

^(A^) = ^{A) , U = const . (4.17) 

Since the global gauge transformations are just rotations in colour space, the wave func- 
tionals have to be colour singlet states satisfying 

L^{A) = (4.18) 

where 

idA'' ^ ' 

is the "orbital" angular momentum in colour space, which is nothing but the colour spin 
of the gauge field. In the next section, we will obtain this constraint in the resolution of 
Gauss' law in the pure Coulomb gauge as a consequence of the zero modes of the Faddeev- 
Popov kernel belonging to the global gauge symmetry. Equation (14.181) implies that the 
wave functional is rotationally invariant, 

^(A) = ^(|A|) (4.20) 

and thus depends only on the gauge invariant modulus |A| of A. 

The global gauge transformations do not exhaust the set of small gauge transformations 
remaining unfixed in the pure Coulomb gauge. An even number of large gauge transforma- 
tions (14. 7p forms a space-dependent small one 

^ / ^ \ '2,7X71 

\/2«(A) = exp [2xK A ■ t) , fe„ = , (4.21) 
which shifts the gauge field by 

A^A^"(^) = (A-2A;„A) . (4.22) 



Note that the gauge transform A^^" still satisfies the pure Coulomb gauge if the original 
configuration A does so. Since the wave functional has to be invariant under small gauge 
transformations, it has to satisfy the condition 



^ A-2A;„A = ^(A) (4.23) 
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and by Eq. (ICTIl 

^(|2A;„-|A||) = ^(|A|). (4.24) 
Restricting A to the first Gribov region 

, . , , 27r , , 

< |A| < A;i = — , (4.25) 

1j 



^(^-|A|) = ^(|A|) (4.26) 



the above condition becomes 

or when expressed in terms of the dimensionless variable "i? = ^-^^ 

^{2n -^) = ^(?9) . (4.27) 

Under large gauge transformations V^(A) (14.71) the wave functional needs only to be invari- 
ant up to a phase 

^ (^A^(A)) = e^"*(A) (4.28) 

and since (V"(A))^ is a small gauge transformation this phase has to be e*" = itl. Using 
Eq. fl4.8p and proceeding as above we find from the effect of the large gauge transformation 
the boundary condition 

^(vr - ^9) = ±*(^9) . (4.29) 

The two signs correspond to two superselection sectors of the theory, which are the discrete 
analog of the 6-vacuum in D = 3+1. In section [6l we will find that the ground state belongs 
to the sector with the plus sign. 



The global gauge symmetry left in the pure Coulomb gauge is used in the diagonal Coulomb 
gauge to diagonalise the (algebra-valued) gauge field 

After implementing this gauge there is still the residual 5*0(2) ~ f/(l) global symmetry 
of rotations around the 3-axis. This abelian symmetry cannot be fixed since the gauge- 
fixed configurations (14.301) are invariant under these rotations. This implies that also the 
wave functional defined on the gauge-fixed manifold automatically respects this symmetry. 
Therefore, the residual global U{1) symmetry can be left out in further considerations. 

A small gauge transformation consisting of two successive displacement transformations 
Vi^s) 114. 8p shifts A^ to A^ — ^. By the identification of±A^ this configuration is equivalent 
to ^ — A^. Thus A^ and ^ — Ac" (or d and 27r — ■(9) represent the same configuration. This 
residual invariance under the small gauge transformations (^(63))^ left by the diagonal 
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Coulomb gauge has to be respected by the wave functional, which therefore has to satisfy 
the boundary condition 

^(27r - = ^(^9) . (4.31) 
This condition (14.311) was already obtained above in the pure Coulomb gauge, see Eq. 
fl4.27p . This is not surprising: Since the diagonal Coulomb gauge contains the pure Coulomb 
gauge the boundary conditions following from the residual gauge invariance in the diago- 
nal Coulomb gauge apply also to the pure Coulomb gauge. In section [6] we will solve the 
Schrodinger equation thereby imposing the boundary conditions fl4.3ip in the diagonal 
Coulomb gauge and the conditions (14.181) and (14.311) in the pure Coulomb gauge. 



5 Resolution of Gauss' law 

As already discussed in section [2], in the Hamiltonian approach in Weyl gauge {Aq = 0) 
Gauss' law ( I2.12p does not follow from the Heisenberg equation of motion and has to be 
imposed as a constraint on the wave functional. In the following, we explicitly resolve Gauss' 
law in both the pure and diagonal Coulomb gauges, thereby paying proper attention to the 
zero modes of the Faddeev-Popov kernel. We will find that the modes = 0, a = 0, ±1 are 
excluded from the Coulomb propagator in the pure Coulomb gauge (although the modes 
n = 0, 0" = ±1 are not zero modes of the Coulomb kernel A„=o,(t±i = 7^ 0!). In the 
diagonal Coulomb gauge only the true zero mode n = a = is excluded from the Coulomb 
propagator. This is in accord with our discussion of the Faddeev-Popov method in section 
m We will first outline the general strategy of resolving Gauss' law and afterwards apply 
it separately to the pure and diagonal Coulomb gauges. 

First note, that even in the pure Coulomb gauge where only constant gauge "fields" are 
left, the momentum operator has space-dependent components. We denote the part of the 
momentum operator conjugate to the modes of the gauge field left after gauge fixing by 
and the remaining part by n||(a;), 

n(x) = nx + n||(x) . (5.1) 

These components are orthogonal to each other in the sense that 

J dxU^U\\{x) = J dxU\\{x)U^_ = (5.2) 

and we refer to them here as the "transversal" and "longitudinal" components of the mo- 
mentum operator, respectively, although this notation is somewhat misleading in the case 
of the diagonal Coulomb gauge (see Appendix [B]). With Eq. ( 15. 2p . the Yang-Mills Hamil- 
tonian (I2.10p becomes 
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As usual, we will solve Gauss' law for the longitudinal part n||(x). Since dll± = (in both 
gauges) we can rewrite Gauss' law (I2.12P as 

D^\x)mx)^{A) = pUx)^{A) , (5.4) 

where 

p^„,(x)=p«(x)+p^" (5.5) 

is the total colour charge, including the external charge p"(a;) and the charge of the gauge 
bosons 

p-^ix) = -A-\x)U', . (5.6) 

Since Gauss' law is a constraint on the wave functional and not an operator identity, one 
can extract from Gauss' law only n||\E'(y4) but cannot obtain flu itself. For this reason, we 
consider the expectation value of the Hamiltonian and perform a partial integration with 
respect to the gauge field to obtain 

2 

= ^ JvaJ dx{U{x)^{A))*U{x)^{A) . (5.7) 

Implementing here the (pure or diagonal) Coulomb gauge by the Faddeev-Popov method, 
splitting the momentum operator into longitudinal and transversal parts 11 = Hn + Il±, 
expressing n||\l/ by Gauss' law and performing a partial integration with respect to the 
gauge-fixed field, the Hamiltonian becomes 

2 

H=^ldxJf^n^JFpU^ + Hc, (5.8) 

where J'pp is the Faddeev-Popov determinant and He is the so-called Coulomb Hamilto- 
nian, defined by 



I)AJpp(A)|-y'c^a:(n||(x)^(A))*n||(a;)^(A) =: Jpp(A)^*(A)ffc^(A) . (5.9) 

Formally, from Eq. (15.41) follows 

mx)^{A) = J dy{x\ {D-'y'\y)pUy)^{A) (5.10) 



and the Coulomb Hamiltonian becomes 

-.2 



Hc = jJ dxdyJ^'ppUx)F^\x,y)JpppUy) , (5.11 
where 

F^\x,y) = {x\{-D-n) = {x 



y) (5.12) 

is the so-called Coulomb kernel. However, the operator D has zero modes, which forbid a 
naive inversion. In Appendix[B], we explicitly solve Eq. (15.41) for nM(x)\E'(y4) and extract He 



23 



for both the pure Coulomb gauge and the diagonal Coulomb gauge, thereby paying proper 
attention to the zero modes. The upshot of these considerations is that the zero modes of 
the Faddeev-Popov kernel, which are a consequence of incomplete gauge fixing, give rise 
to additional constraints on the wave functionals. These constraints basically arise from 
the projection of Gauss' law onto the zero modes of the Faddeev-Popov kernel. In the pure 
Coulomb gauge these constraints read (IB. 311) 

Q^vP = , (5.13) 

where 

L L 

Ql, = J dxpU^) = jdx {p\x) + pl{x)) = g» + Ql (5.14) 



is the total colour charge. In this gauge the transverse momentum operator flA.TP reads 



^ d , , 

ni = T-nr 5.15 

^ Lid A'' ^ ' 

and the dynamical charge of the gauge bosons (EEj) becomes (up to a sign) the colour 
angular momentum operator (more precisely the colour spin) of the gauge field, (14.191) . 
The residual constraint ( IB. 131 ) from Gauss' law becomes 

= Q'^^ , (5.16) 

where is the external charge. In the absence of external colour charges = this 
constraint simplifies to 

= , (5.17) 

i.e. the physical wave functionals do not depend on the angle degrees of freedom A(6', 0), 
which, in fact, are unphysical since they represent the residual global colour gauge degrees 
of freedom, which are not fixed by the pure Coulomb gauge condition. Thus, the physical 
vacuum wave functionals depend only on the "radial" coordinate |A| which is the physical 
degree of freedom of the gauge field. The constraint ( 15.17^ was already found in the previous 
section, see Eq. ( I4.18p , and refiects the invariance of the wave functional under global gauge 
transformations. 

In the diagonal Coulomb gauge, where 

A" = S^'^A^ , ni = S^'^Ul , (5.18) 

the dynamical charge of the gauge bosons (15.61) vanishes and the residual constraint 
from Gauss' law implies the vanishing of the Cartan component of the external charge in 
the physical state (see Eq. ( IB. 321 )) 

Q^^ = 0. (5.19) 

In absence of external charges = 0, in this gauge there is no residual constraint on the 
wave functional from Gauss' law. 
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The diagonal Coulomb gauge rotates the constant gauge mode in the 3-direction, i.e. = 
±|A| and with the restriction to the fundamental modular region (14.151) . in this gauge 
equals the modulus of A. Thus, in both gauges the physical wave functional depends only on 
the modulus of the constant gauge mode A"-, which is the only physical degree of freedom. 
Both gauges leave a residual global gauge invariance: Global SU (2) symmetry in the case 
of the pure Coulomb gauge and global U{1) symmetry in the case of the diagonal Coulomb 
gauge. By Noether's theorem, these global symmetries imply the existence of conserved 
charges: Qm^''^'^ in the case of the pure Coulomb gauge and in the diagonal Coulomb 
gauge. The residual constraints on the wave functionals obtained above from Gauss' law are 
nothing but the quantum version of Noether's theorem for these global colour symmetries. 
Also in 3 + 1 dimensions the pure Coulomb gauge fixing still leaves invariance with respect 
to global colour gauge transformations and by Noether's theorem the total colour charge 
has to be conserved [351. 



After resolution of Gauss' law (see Appendix [B|) one finds the following gauge-fixed Hamil- 
tonian in the pure Coulomb gauge 

^ d ^ d , , 

where Jp is defined by Eq. (13.32^ and the Coulomb Hamiltonian (cf. Eq. f l5.1ip ) is given 
by 

2 

Hc = yI d^dy p\x)F^\x, y)p\y) (5.21) 

with Coulomb kernel 

F-\x,y) = F'^'IA = AX](x,y) = E(^l^) T.i^\UW)\-'A^\U'^\b){n\y) , (5.22) 

from which all zero modes = of the Faddeev-Popov operator (I3.24p are excluded, 
although n = 0,a = ±1 are not zero modes of the operator (— -D^) in Eq. (15. 12^ ! Note 
also that the dynamical charge pg (15. 6p . although being here non-zero, has dropped out 
from the Coulomb Hamiltonian (I5.2ip . This is a special feature of 1 + 1 dimensions (see 
Appendix [B]) and is a consequence of Pg being space-independent in the Coulomb gauge. 
Similarly, the Faddeev-Popov determinant also drops out from the Coulomb Hamiltonian. 

The first term in Eq. ( I5.20p arises from the "transversal" momentum operators cor- 
responding to the physical mode {A°- = const). This term has the form of a Laplacian 
in a curved space with the Faddeev-Popov determinant acting as the determinant of the 
metric. The second term of Eq. ( 15.20p arises from the "longitudinal" (here x-dependent) 
part n||(a;) of the momentum operator. This term gives the static potential of external 
static colour charges, and it is considered an advantage of the pure Coulomb gauge that 
this term is explicitly isolated. Note, the Hamiltonian in the pure Coulomb gauge (15.201) is 
still invariant under global colour rotations, which are not fixed in this gauge. 
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Resolving Gauss' law in the diagonal Coulomb gauge (I2.20p (which does fix the global 
colour rotations) yields the gauge-fixed Hamiltonian 



where J7d is defined in Eq. f l3.20p . Here the Coulomb Hamiltonian He is still given by Eq. 
(I5.2ip . however, with the Coulomb kernel F'^'^[A]{x,y) (15 .22^ replaced by 

F^'ix^y) = F'^'IA = \A\T,]{x,y) =J2'{^\n){a\a)X;,l{a\b){n\y) , (5.24) 



where the prime indicates that the mode n = cr = is excluded, which is the only zero 
mode of the Faddeev-Popov operator in this gauge. 

Note, since the pure Coulomb gauge field A is related to the field in the diagonal Coulomb 
gauge by a global gauge transformation, see Eq. (13.251) . one would expect that in view of 
Eq. (13.26p . the Coulomb kernels (15.121) in these two gauges are related by 

F^'iA = UA'T,U^{x,y) = UacF''[A'n]ix,y)Uj, . (5.25) 

This is almost the case (cf. Eqs. (I5.22p and (I5.24p ) except for the additional zero modes 
n = 0, 0" = ±1 to be excluded from the kernel (I5.22p in the pure Coulomb gauge. The zero 
modes n = 0, cr = 0, ±1 are precisely given by the constant gauge transformation (I3.29P 
= {a\U\a){{x\n = 0) = const) not fixed in the pure Coulomb gauge. 

By fixing the pure Coulomb gauge, one switches from Cartesian to curvilinear coordi- 
nates and accordingly the gauge-fixed Hamiltonian acquires the form of the Hamiltonian 
in curved space The gauge-fixed Hamiltonian is, of course, no longer gauge invariant. 
In particular, the Hamiltonian (15.23P is not even invariant under global colour rotation 
since the diagonal Coulomb gauge (12.201) fixes also the global gauge transformations. Nev- 
ertheless, the diagonal Coulomb gauge still leaves invariance under global abelian gauge 
transformations (colour rotations around the 3-axis). 

In the absence of external charges {p"'{x) = 0), the Coulomb term He (I5.2ip obviously 
vanishes in both gauges, however, for different reasons, see appendix[Bl In the pure Coulomb 
gauge it vanishes because all constant modes n = 0, a = 0, ±1 are excluded from the 
Coulomb kernel F"-^{x,y), while in the diagonal Coulomb gauge it vanishes because the 
physical modes live all in the Cartan algebra, resulting in a vanishing colour charge (15. 6p 
of the gauge bosons. 
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6 The physical state space 



The spectrum of the Yang-Mills Hamiltonian in 1 + 1 dimensions can be obtained in a 
gauge invariant way |36| . We choose here the diagonal Coulomb gauge, with all unphysical 
degrees of freedom eliminated and the physical degree of freedom within the fundamental 
modular region < -d < tt, see section [H Thus, one can regain the gauge invariant spectrum 
and simultaneously find the physical state space. From the vacuum wave functional "^[A] 
in the diagonal Coulomb gauge, the one in the pure Coulomb gauge can be derived which 
will be very useful in the subsequent sections. 



6.1 Diagonal Coulomb gauge 



In the diagonal Coulomb gauge (I2.20p . the Yang-Mills Schrodinger equation of the 1 + 1 
dimensional theory 

F^'fc = Ek^k (6.1) 

can be solved exactly [28]. In the absence of external colour charges the Yang-Mills Hamil- 
tonian (I5.23P reads in the compact variable = ^-^^ 

1 (gL\^ 1 d ^ , d 



To solve the Schrodinger equation, we introduce the "radial" wave functional 0(^9) by 

'^{^) = ^=m- (6.3) 



This eliminates the Faddeev-Popov determinant Jd{,'&) = sin^ in the scalar product 

TT TT 

= / d^JD{^)^l{^)^2{^) = J d^(j)l{^)M^) (6.4) 



and reduces the Schrodinger equation (16.11) to 

The boundary condition (14.311) on the total wave functionals ipi'^) requires the radial wave 
functional 0(^9) to satisfy 

0(27r -^) = -0(^9) . (6.6) 



27 



Thus, the solutions to Eq. (16.5P read 

^^(^) = ^ sm{k^) , kEN. (6.7) 

These are normalised with respect to the scalar product (16 ■4p . The corresponding energy 
eigenvalues are given by 

= '-^{k' - 1) = ^ J(J + 1) = ^ 1(1 + 2) . (6.8) 

Here, we have defined 

k-1 1 , , 

J = ^ = 0, 2^ 1, ••• (6.9) 

to identify the spectrum (16.81) as a rigid rotor in colour space where the integer and half- 
integer j correspond to the two superselection sectors defined by the boundary condition 
(I4.29p . Alternatively, one can use the definition 

l = k-l = 2j (6.10) 

to recognise in Eq. f l6.8p the energy eigenvalues of a point particle with mass 4:/{g'^L) and 
angular momentum / on a unit sphere S3 in D = 4, which is the group manifold of SU{2). 
(In fact, H (16. 2p is (up to the constant factor) the polar angle part of the Laplacian on 
5*3.) Either way, the eigenfunctions 

/2'sin((27 + 1)?9) [2 , , , 

= \ = \-Xjm 6.11 

V smw V ^ 

are the characters of SU (2) (fi - arbitrary unit vector) 

X,(/9)=E0'^|e^^"1j'^) (6-12) 

TO 

and the eigenvalues (16. 8p are seen to diverge in the thermodynamic limit L ^ 00 except 
for j = 0. Therefore, all states are frozen, except the one with j = 0, which has vanishing 
energy {Ek=i = 0) and the vacuum wave functional 

^(^) = \/f- (6.13) 

The vacuum wave functional \l/ has no dependence on the gauge "field" t9 = -^^^ and 
describes a stochastically distributed weight of gauge field configurations in the expectation 
values of the diagonal Coulomb gauge, 

TV rt 

(C(^))diago„al gauge = / Jn{^)0{^)\^ = -J d^Sm^^O{^) . (6.14) 
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It is claimed that in 3 + 1 dimensions the gauge field configurations that dominate the 
infrared physics are located on the common boundary of the Gribov region and the fun- 
damental modular region [37]. In 1 + 1 dimension, this is obviously not realized. The 
Faddeev-Popov determinant Jr, = sin^ actually suppresses the contributions of "i? = vr 
which is the common boundary of Gribov and fundamental modular regions. This is due 
to the fact that in the fully gauge-fixed theory, i.e. in the diagonal Coulomb gauge, the 
configuration space is merely one-dimensional and entropy does not favour boundary con- 
tributions (although expected to do so in D = 3 + 1). 



6.2 Pure Coulomb gauge 



We are interested here in the vacuum wave functional in the pure Coulomb gauge. Although 
the pure Coulomb gauge Hamiltonian (I5.20p apparently has a much more complicated 
structure than the Hamiltonian (16. 2p in the diagonal Coulomb gauge both expressions 
are equivalent, due to gauge invariance. In the following, we will explicitly reduce the 
Hamiltonian in the pure Coulomb gauge to the diagonal Coulomb gauge Hamiltonian (16.21) . 
This will provide us with the explicit form of the wave functional in the pure Coulomb 
gauge, which is needed for subsequent considerations. 

In the absence of external charges, the Yang-Mills Hamiltonian (I5.20p is proportional to 
the Laplacian Ac in the space of gauge orbits projected on the hyperplane defined by the 
pure Coulomb gauge condition 

H = -^-^Ac, Ac := j-V ■ JpV , V = e„^. (6.15) 

Using spherical coordinates of the gauge "field", introduced already in Eq. (12.251) . 



A = AA(0,0), A=\A\ = VA^, (6.16) 
and expressing the V operators in Eq. (16.15P by the standard form 

where L is the colour angular momentum operator (14.191) . the Laplacian Ac reads 

_ 1 d d 

Here, we have used the explicit forms of the Faddeev-Popov determinants J^p in the pure 
Coulomb gauge (13.32^ and J^c in the diagonal Coulomb gauge (I3.20p . In the physical space 
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of colour singlet states where (see Eq. (I5.17P ) 



L2|^) = 0, (6.19) 

the last term in Eq. ( 16.18^ becomes irrelevant and the Hamiltonian H in Eq. ( IG.lSp reduces 
precisely to the one in the diagonal Coulomb gauge f l6.2p . Its eigenfunctions in the colour 
singlet Hilbert space (I6.19P are therefore the same as in the diagonal Coulomb gauge. 
However, in the pure Coulomb gauge the global gauge degrees of freedom 6, defining the 
orientation A(5, 0) of the colour vector A remain as coordinates, which enter the definition 
of the scalar product in the Hilbert space of the wave functionals, 

(^i|^2) = ^JdQj d^^^Jp{^)^l{^,e,<P)<i/2{^,0,<l)) ■ (6.20) 

52 

Here, dfl{9, 0) denotes the usual integration measure on 5*2 and the Jacobian 'd'^ comes from 
the transformation to spherical colour coordinates. Due to the factor in the definition 
of the scalar product (16.20p . the normalisation of the vacuum wave functional in the pure 
Coulomb gauge is identical to the one in the diagonal Coulomb gauge in Eq. fl6.13p . 



^(A) = ^- . (6.21) 
The pure Coulomb gauge expectation values 

{0[A]) = JvAJp[A]^*[A]0[A]^[A] 

TT 

d^^^Jp{^)^j c/fiC[^9,^,0] 

52 

j d^ sin^ ^— j dn 0['d, 9, 0] (6.22) 



TT 
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average the gauge degrees of freedom 9, 0. The gauge invariant variable is integrated with 
the same weight as in the diagonal Coulomb gauge expectation value f l6.14p . For operators 
0{'d) in the diagonal Coulomb gauge, the expectation value (16.22^ gives the same result 
as the one in the diagonal Coulomb gauge (16.141) . Hence, we may use the definition (I6.22P 
exclusively. 

In subsequent sections we will use the exact wave functionals to calculate various propa- 
gators and vertices, the colour Coulomb potential and derive their Dyson-Schwinger equa- 
tions. 
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7 Exact Propagators and Vertices 



In D = 3 + 1, the study of Landau gauge Dyson-Schwinger equations (DSEs) has recently 
become quite popular (for a recent review, see Ref. [38] and references therein). The co- 
variant Landau gauge is technically convenient and the structure of the DSEs is similar to 
the pure Coulomb gauge. In the pure Coulomb gauge, the variational approach to D = 3 + 1 
Yang-Mills theory also results in a set of DSEs for the propagators and vertices |H|. In any 
case, a truncation of the non-linearly coupled DSEs is unavoidable. It is difficult to assess 
the validity of this approximation. In 1 + 1 dimension, on the other hand, we can calculate 
the exact Green functions and compare them to the solution of approximated DSEs. In this 
section, we will use the exact vacuum state (I6.2ip to calculate the propagators and vertices 
in the pure Coulomb gauge. We begin with the gluon propagator. After the computation 
of the ghost and Coulomb propagators, the ghost-gluon vertex is calculated in the pure 
Coulomb gauge. Finally, we discuss the propagators in the diagonal Coulomb gauge. 



7.1 Gluon propagator in pure Coulomb gauge 

In 1 + 1 dimensions, the pure Coulomb gauge fields are spatially independent. The gluon 
propagator is therefore a constant matrix, defined as the expectation value of two field 
operators. 

Expressing the colour vector A"- = j^'d A."- {6 , (p) in spherical coordinates, the pure Coulomb 
gauge expectation value (16.221) for the gluon propagator (|7.1I) yields 

2 ^ 

D''^ = -(-\ I dd sin^ d — I dVL A%e, <f))A\e, <f)) . (7.2) 



By symmetry, the angular integration yields 

j-JdQA%e,<p)A\d,<p)=^-6''' (7.3) 

and thus for Eq. ([72]) 

j^ab ^ ^ab}_ (^^2 _ 2^ ^. ^ab ^ ^ 3 1. . (7.4) 

The pure Coulomb gauge gluon propagator has only diagonal components which all have the 
same value D^. In the thermodynamic limit, L — oo, the gluon propagator is identically 
zero, in agreement with the fact that the theory becomes trivial for L oc. 
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1.2 Ghost propagator in pure Coulomb gauge 



The ghost propagator occurs in the Dyson-Schwinger equations as a consequence of the 
projection on the hypersurface of Coulomb gauge (or Landau gauge). In the variational 
approach [8], this propagator is merely an auxiliary object to facilitate the computation 
of the energy density. It is defined as the expectation value of the inverse Faddeev-Popov 
kernel. The Faddeev-Popov kernel Ai of the pure Coulomb gauge is given in coordinate 
space by Eq. (]3.24p . In the momentum space, we have 

L 

Mt = J dxM''\x,y) e-""-' = - tA'^'kn =: (G;^)"' . (7.5) 





We refer to the inverse M.^^ as the "ghost kernel" denoted by Gn- It is customary to 
consider the Cartesian colour components of the matrix- valued ghost propagator. Let us 
therefore invert A^„ in the Cartesian basis. Using the SU{2) identity 

{AAy^ = A^'A^ - S^^A^ , (7.6) 

one can verify that 

G? = ^^U°'-^ + 4 (7.7) 

is indeed the inverse of ^A'^ (17. 5p . The expectation value ()6.22l) then defines the ghost 
propagator 

[Gt) = ^ (7.8) 

as well as the ghost form factor d^, which measures the deviation of {G'^) from the tree- 
level behaviour 6°"^ /k'^. The angular averages within the pure Coulomb gauge expectation 
value (17. Sp can be taken with Eq. (17.31) and the identity 

I dnA''^ = Q. (7.9) 



A-n 

One thus finds a colour diagonal ghost form factor, 

/ 1 / 1 A2\ V 

lab xab / i i \ \ . ):ab 



Its diagonal elements dn can be rewritten as ( ^ = — 



dn = l + -(^^, ^) = 1 + T^ d'&sm^^^— -. (7.11) 



3 \ (nn)'^ — "d"^ Sn J (nn)'^ — 1?^ 
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Fig. 3. The exact ghost form factor, given by Eq. (|7.1ip . 

The above integral may be expressed by integral sine functions. The allovv^ed modes n G 
Z \ {0} exclude the zero mode n = of the Faddeev-Popov operator. In the ultraviolet 
limit, the ghost form factor (IT.lip approaches tree-level, 

lim d„ = 1 , (7.12) 

n— >oo 

since the D = 1 + 1 theory is super-renormalisable and there are no anomalous dimensions. 
In the infrared, the ghost form factor is enhanced, as can be seen in Fig. [3l This enhance- 
ment is also found in D = 3 + 1 dimensions and is understood to come from near-zero 
eigenvalues of the Faddeev-Popov kernel in the vicinity of the Gribov horizon. 

In the 3 + 1 dimensional continuum theory, the inverse of the ghost form factor, 
represents the generalised dielectric function e of the Yang-Mills vacuum [39j. Fig. [3] shows 
that also in D = 1 + 1 the Yang-Mills vacuum behaves like a dia-electric medium (e < 1) 
in the infrared and becomes the ordinary vacuum (e = 1) in the ultraviolet. 



7.3 Coulomb propagator in the pure Coulomb gauge 

The pure Coulomb gauge Hamiltonian (I5.20p comprises the so-called Coulomb term He, 
see Eq. f l5.2ip . which accounts for the interaction energy between colour charges. By the 
calculation of the expectation value (-ffc), the quark potential can be found, see section 
m The "Coulomb kernel" F, which mediates this interaction is defined in Eq. (|5.22l) and 
reads in momentum space 

f:' = Gi^kic^:: . (7.13) 

Using the explicit form (I7.7P of the ghost kernel and the identity (17.61) . the Coulomb 
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Fig. 4. Form factor of the Coulomb potential as given by Eq. (I7.17p . 
kernel F^'' is cast into the form 




(7.14) 




The expectation value (16 .22^ of the operator F^* in Eq. (17.141) defines the Coulomb prop- 
agator 

{f:") = ^ (7.15) 

and the form factor v^hich measures the deviation of {F^'') from tree-level, S"'^/k'^ (being 
the abelian case). To evaluate the expectation value (17.15^ . let us first integrate the gauge 
degrees of freedom of F^^{A). Using the identities (17. 3p and (17. 9p . one finds from F^^ in 
Eq. (jmi) 



— / dQFf 
47r 



1 / ^ 



i + f-V + ^f-r f-r-3 



2 I .q2 



(7.16) 



The form factor (17.150 of the Coulomb propagator is therefore strictly diagonal and its 
diagonal elements yield 

=: 5«Vn , 0„ = - I 1 + 47^2 / (I'd sin^ ^ {Trnf + ^^ \ ^ ^ 

In the ultraviolet limit, the form factor 0„ approaches tree-level. As shown in Fig. [H 0„ 
is infrared enhanced. In the D = 3 + 1 theory, the infrared enhancement of the Coulomb 
propagator is expected to come from the restriction of the configuration space to the Gribov 



34 



g I I I I I I I I I I 

1 23456789 10 
n 

Fig. 5. Coulomb form factor /„, as defined in Eq. (17.181) . 

region, as claimed by Gribov in his seminal paper [40], and to lead to a confining quark 
potential. This will be discussed further in section [3 

An issue in D = 3 + 1 dimensions is whether the expectation value of (I7.13P can be fac- 
torised, i.e. can the connected part be neglected? In order to answer this question, a further 
form factor was introduced in Ref. [Sj. It measures the deviation from the factorisation, 

(G„ klGn) =■■ {Gn) kl U (Gn) , ^ /n = ^ • (7.18) 

and can be expressed by the ratio of the form factor 0„ of the Coulomb propagator and the 
ghost form factor squared. Following Ref. [8], we refer to /„ as the "Coulomb form factor". 
In 1 + 1 dimensions, where the exact solutions for dn and 0„ are available, we can calculate 
fn exactly. In Fig. [5l the result for is depicted. It shows an infrared enhancement. A 
further discussion in the context of the D = 3 + 1 theory will follow in section [121 



7.4 Ghost- gluon vertex in pure Coulomb gauge 

In the variational approach in the pure Coulomb gauge ^ and the Dyson-Schwinger ap- 
proach in Landau gauge, the ghost-gluon vertex is of particular interest. In these ap- 
proaches the proper ghost-gluon vertex is usually replaced by the bare one with the argu- 
ment that this vertex is not renormalised [H]. In fact, recent lattice calculations performed 
in D = 1 + 1 Landau gauge provide little evidence for a dressing of this vertex [42]. However, 
although larger lattices are available in D = 1 + 1, there are significant statistical errors. 
The lattice results are also plagued by the existence of Gribov copies. The 1 + 1 dimensional 
continuum theory, on the other hand, has full control of the Gribov problem, see section [H 
In the following we will calculate the proper (one-particle irreducible) ghost-gluon vertex 
in the pure Coulomb gauge. 
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The bare ghost-gluon vertex is defined by [8] 

T''^{x, x') = -'^^^j^ = f'^dJix, x') , (7.19) 
where G^^ = M. in the Faddeev-Popov kernel ( 13.24p . Fourier expansion 

r°''^(x, x') = \y. e^'^-^^-^'^rO'" (7.20) 
yields the momentum space representation 

r°''^ = iknf" . (7.21) 

In the D = 1 + 1 Coulomb gauge the ghost-gluon vertex depends only on a single momentum 
for there is only the zero momentum mode of the gauge field. The proper gho st-g luon vertex 
is defined via the expectation value for the connected ghost-gluon vertesLfJ 

{A^G':) = D^^' (c') ixif {Gi^) = da^ i^r ^ (7.22) 

which is the proper vertex with propagators attached on its legs. We have used the fact 
that the exact ghost and gluon propagators in the pure Coulomb gauge are colour diagonal. 
Using the explicit form (17.70 of the ghost kernel G^, we can write the expectation value 
(17:221) as 

1 / 1 / 4* 4^^ 4^^=^ \ 



A-G^:) = - I^A-r - A^^ + ^A-— j j . (7.23) 



k 



In the pure Coulomb gauge expectation value ( 16.22p . the angular integration renders the 
first two terms in Eq. ( 17.23^ zero, while the last one yields 

1 / A"- A'^ \ 

A''G':) = -ikniT^'''^ \ 



1 1/792 

-*- f-pO,d\bc ^ i-ad I 

A;2 ^ " ^ 3 \ (7^)2 - 

k 



A^nT^{dn-l) (7.24) 



where in the last line we have used the expression (I7.1ip for the ghost form factor dn. It is 
helpful to define the form factor 7„ of the ghost-gluon vertex by 

K=:7nTy. (7.25) 



^ For the wave functional chosen in Ref. [8], the definition (|7.22|) coincides with the one used 
there, (GT^G) =: (G)r(G). 
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Fig. 6. Form factor 7„ of the proper ghost-gluon vertex, as given by Eq. (|7.26l) . 
This form factor can nov^^ be expressed by Eq. (17.241) using the definition (I7.22P of 

k'^{dn 1) 



In 



(7.26) 



Let us check the ultraviolet limit, n ^ oo. From Eqs. f l7.1ip . f l7.2p and (I7.3p . one can see 
that 

16 } 

lim klidn - 1) = — / sin^ ^ = 2D A (7.27) 



where Da is the gluon propagator defined in Eq. f l7.4p . It then follows immediately from 
Eqs. f l7.26p and (I7.12p that the ghost-gluon vertex approaches tree-level asymptotically, 



lim 7„ 



1 . 



(7.28) 



In Fig. El one can see that the form factor 7„, given by Eq. ( 17.26^ . hardly deviates from 
tree-level for the entire momentum range. Deviations are in the range of 5%. This strongly 
supports the popular approximation of truncating Dyson-Schwinger equations by choosing 
a tree-level ghost-gluon vertex. Further discussion will follow in section [TTl 



7.5 Propagators in diagonal Coulomb gauge 



For a comparison to the pure Coulomb gauge, which does not completely eliminate the 
gauge degrees of freedom, we here calculate the propagators in the diagonal Coulomb 
gauge with restriction to the fundamental modular region. The gluon propagator Z)"^* in 
the diagonal Coulomb gauge is most easily found by letting A°- ji!} S"-^ before taking the 
expectation value (17.10 . We can then express D"^'^ in the diagonal Coulomb gauge by the 
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pure Coulomb gauge result (I7.4p for Da, 



D'' = 1[y) Jd^^'sm'^ = 3DA, D'' = D'^ = (7.29) 



The colour trace of the gluon propagator is evidently the same in both gauges. From the 
observation that the diagonal Coulomb gauge can be reached from the pure Coulomb gauge 
by the unitary transformation (see appendix [A|) 

A^UAU^ (7.30) 

of the matrix- valued variable A = A"'T°', the invariance of the colour trace of the gluon 
propagator is seen to be an immediate consequence. 



The Faddeev-Popov kernel A^„ in the diagonal Coulomb gauge is substantially different 
from the one in the pure Coulomb gauge, see Eq. f l3.9h (cf. Eq. (13.24^ ). It is convenient to 
expand Ain in the spherical basis ( 1A.13I ) where Ain is diagonal, 

(a'lMnW) = An,Ja'a =■■ (7.31) 



and the eigenvalues A„ o- are given by Eq. (13.161) . We use the vacuum expectation value of 
the ghost kernel in the diagonal Coulomb gauge (I7.3ip to define the ghost form factor 
components (iJ^jZl 

a = ±l 

{GD =: ]: . (7.32) 

The vacuum expectation values (|6.22l) then yield for 7^ 

<=^' = (^-^), <=° = 1- (7.33) 

Taking the colour trace J2a dn of the ghost form factor in the diagonal Coulomb gauge gives 
the same result as summing the diagonal elements d'^°' of the ghost form factor in the pure 
Coulomb gauge (see Eq. (17. lip ). 

E< = 1 + 2 (737X^1) = EC • (7.34) 



^ For the components a = ±1, we related the form factor to the propagator differently from 
the cr = component, in order for all d^ to have the same dimension and complex phase. This 
circumstance which does not occur in the pure Coulomb gauge can here be traced back to the 
structure of the gauge condition (|2.20p of the diagonal Coulomb gauge. 
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While the abelian component a = of the diagonal Coulomb gauge ghost form factor is at 
tree-level, see Eq. (17.330 . the other diagonal components are larger than the pure Coulomb 
gauge result (17.111) for dn, such that the colour trace is invariant. 

The same scenario occurs for the Coulomb propagator (Fn). With the the definition ( I5.24p 
of Fn and Xn,a given by Eq. ( I3.15p . we have with F„|cr) = F^\a) 

and find for the form factor 0^ of the Coulomb propagator for n 7^ 

Noting that the Coulomb kernel F^ in the diagonal Coulomb gauge actually follows from 
a rotation in colour space from the one in the pure Coulomb gauge, see Eq. ( I5.25p . the 
invariance of the colour trace is clear. 

The above results exclude the n = modes which are the zero modes of the Faddeev- 
Popov operator in the pure Coulomb gauge. In the diagonal Coulomb gauge, however, the 
constant n = modes with a = ±1 are allowed. These give the results 

(g;z„") = (i) K 0.39 iaL (7.37) 

(f„"±') = ^(l)«0.23L' (7.38) 

It turns out that these values are of no importance for the considerations below and are 
just given here for completeness. 



8 The static quark potential 



The gauge invariant potential energy of a quark-antiquark pair in = 1 + 1 is well-known 
from the calculation of the temporal Wilson loop. A linearly rising potential emerges and 
the corresponding string tension aj shows strict Casimir scaling [43] . 

= + (8.1) 

In the fundamental representation we have with j = \ the string tension 

a:=a._ = lg\ (8.2) 
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In the Hamiltonian approach, the static colour potential can be obtained by taking the 
expectation value of the Hamiltonian if the charge distribution p^{x') is chosen to be a 
pair of opposite point charges of strength and — g'^, localised at the positions x and 
respectively, 

p»(a;') = {5{x' -x)- 5{x' - y)) . (8.3) 

This calculation is shown here to also give rise to a linear potential which is, however, a 
gauge dependent quantity. Therefore, the string tension ac defined by this linear potential 
gives only an upper bound [37f44] on the gauge invariant string tension a in Eq. (18.21) . 
Nevertheless, the Coulomb potential is calculated here for comparison to the potential 
calculated in the D = 3 + 1 theory. 



8.1 Pure Coulomb gauge 



We first fix only the pure Coulomb gauge. For the colour charge density (|8.3p , the Coulomb 
Hamiltonian ()5.21l) reduces to 
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He = Y^^q' [F'^ix, x) + F''\y, y) - F''\x, y) - F^\y, x)\ , (8.4) 

with the Coulomb kernel F"''(x, y) f l5.22p in the pure Coulomb gauge. The first two terms of 
Eq. (18. 4p represent the self-energy of the static point charges. Only the charges belonging 
to the generators of the Cartan subalgebra can be specified and the expectation value of 
(18. 4p in the Yang-Mills vacuum state for abelian (Cartan) unit charges defines the static 
quark potential. With = 5"^ for SU{2)^ the static quark potential Vcij-) with r = \x — y\ 
becomes 

Vc{r) = {Hc\^^s^, = g' {F'%x,x) - F''ix,y)) (8.5) 

To obtain the above expectation value, we first take the angular average of the momentum 
space Coulomb kernel F^*, see Eq. (I7.16p . 



^Jdn {F''ix,x)-F''{x,y) 

= -y — [ dQF^' (i 



e 

t - i^-^y f 1 + 2i.nr + (1 - cos n)) . (8.6) 
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With 1 — cos(2a;) = 2sin^(a), the Fourier transformation yields 



Air 
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1 L ^ / 1 ^ d ^ \ ■ 2/.r \ 



It + " " (2 ^ ^) - - ^) (2 i ^))) (8-7) 



3 V2 

Here, we have used the formulae in Ref. [45] to obtain for the above sums 

- sin^ (f n) ^ sin(^^) sin(ft^) ^ vr^ / r _ x r x 2 
n2 - (^9/7r)2 2 ^9sin?9 ' 2 VL 

The (elementary) "i^-integration of the expression (18.71) as defined in the expectation value 
(16.22P yields the static quark potential ( IS.Sp 

/ N f . 1 ^ 2r - L sin(27ry)\ / ,^/r\\ , , 

fcM = |(L + 5.-5^ + i — ^)=^.(i + oy) . (8.9) 

In the thermodynamic limit, the string tension ac is defined by 

lim Vc{r) = acr , ^ ac = ^ > cr (8.10) 

L— >oo J, 

and is found to be larger than the gauge invariant string tension a (18.21) . as expected. 



8.2 Diagonal Coulomb gauge 



Let us now apply the global gauge transformation from the pure to the diagonal Coulomb 
gauge and recalculate the string tension of the potential. The 33-component of the Coulomb 
kernel F^'' can be found by setting A" = j^i!}6'^^ in Eq. (17.141 ). 

= ^ . (8.11) 



It apparently mediates only the tree-level (abelian) part of the Coulomb interaction. Since 
F^^ is thus field independent, the vacuum expectation value needs not to be taken and we 
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Coulomb gauge 
diagonal gauge 




Fig. 7. Coulomb potential Vc{r) for the pure and the diagonal Coulomb gauge. 



directly arrive at the Coulomb potential (I8.5P by Fourier transformation 
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Voir) = / y E (l - ^-''"') = / 4 E 



n=l 



2 



r 1 + C 



(8.12) 



The string tension of the diagonal Coulomb gauge is obviously identical to the one in the 
pure Coulomb gauge., see Eq. (18.101) . It can actually be seen by taking the limit L ^ oo 
in the pure Coulomb gauge expression (18. 7p that there Coulomb kernel also turns out field 
independent (i.e. ^ independent) and coincides with the Fourier transform of the diagonal 
Coulomb gauge kernel ( 18. lip . 

The potential Vc{r) is shown in Fig. [7] for the result in the pure Coulomb gauge (18.91) and 
in the diagonal Coulomb gauge (18.121) . It is clear that in either gauge the function Vc{r) 
must be symmetric about the axis r = ^, since a separation of charges by r is identified 
with a separation by L — r on the spatial manifold S^. For <^ 1, the potentials Vc{r) 
are seen in Fig. [7] to have the same slope in both gauges, i.e. the string tensions have the 
same value ac given by Eq. (18.100 . 



Strictly speaking, the sum of momentum modes n must be changed to which excludes the 
zero mode n = cj = of the Faddeev-Popov operator in the diagonal Coulomb gauge, but includes 
the modes n = and o" = ±1, see section [3l The result for these zero modes is given in Eq. (I7.38p . 
However, in the subtraction of eigenenergies in the potential Vc{r), see Eq. (|8.5p . all modes with 
n = cancel and it makes no difference whether Y^^j^q or is used. 
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9 Dyson-Schwinger equations 



In this paper, it is intended to test the approximations made in the study of Dyson- 
Schwinger equations (DSEs) in higher dimensions by considering the D = 1 + 1 case. The 
exact DSEs for the ghost and gluon propagators are usually derived from the partition 
function of Yang-Mills theory. In the present case, we can simply use the definition (I6.22P 
of vacuum expectation values to come by this set of equations. We restrict ourselves to 
the exact ground state of the pure Coulomb gauge, leaving aside the global rotation to the 
diagonal Coulomb gauge. It will be shown how Gribov copies affect the Dyson-Schwinger 
equations and their solution. 

Let us start with the derivation of the DSE for the gluon propagator. It follows directly 
from the expectation value f l6.22p that {^{A) = const) 

° = / ^•^il^-^"'^' MAW-e^'^' = ((^ + ^°) (9-1) 

holds, since the integral of the total derivative in Eq. (19. ip is proportional to the Faddeev- 
Popov determinant J7p evaluated at the first Gribov horizon dQi, where it vanishes. Ap- 
plying a derivative S/6j'' and setting the sources j'^ to zero gives 



= ^A' -jjf- + y = - Tr G r + b"'' . (9.2) 

The trace "Tr" in Eq. (19. 2p sums up the diagonal elements in colour space as well as all 
modes kn with n = ±1,±2 . . . , excluding the zero mode n = of the Faddeev-Popov 
determinant 

Jp = expTrlnG"^ (9.3) 

in the pure Coulomb gauge. We recognise in Eq. ( 19. 2p the connected ghost-gluon vertex 
{A'^G). With its decomposition ( I7.22p into the proper vertex and the attached propa- 
gators in momentum space, Eq. (19. 2p can be written after contracting with 5"^ as 

(9.4) 




This Dyson-Schwinger equation holds for any SU{Nc) but we eventually set Nc = 2 for 
comparison with the exact SU{2) results in chapter [71 Using the definition ( I7.25p of the 
form factor 7„ for the ghost-gluon vertex and tr T^T*^ = —Nc{N^ — 1), the DSE for the 
gluon propagator (19. 4p can be written more concisely. 

Da' = E 7n § . (9.5) 
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Fig. 8. DSE for the gluon propagator. Curly lines represent the gluon propagator, dashed lines 
the ghost propagator. The empty blob represents a proper ghost-gluon vertex, the dot a tree-level 
ghost-gluon vertex. 



-1 -1 



Fig. 9. DSE for the ghost propagator. The dotted line is the tree-level propagator. 
Diagrammatically, the inverse gluon propagator is given by a ghost loop as shown in Fig. 



[Hi ^^1 Inserting into the right-hand side of Eq. (|9.5I) the exact expressions for the ghost 
form factor dn (IT.lip and the ghost-gluon vertex 7^ (17.261) . one can explicitly show that 
the ghost loop (r.h.s. of Eq. ( 19. Sp ) equals the expression for Dj^ obtained in Eq. (17. 4p . In 
Ref. pj, the ghost loop is referred to as the "curvature" since it incorporates the curvature 
of the space of gauge-fixed variables. It is found that the curvature governs the infrared 
behaviour of the gluon propagator such that the exact DSE of 1 + 1 dimensions in Fig. [8] 
is the infrared limit of the corresponding DSE in 3 + 1 dimensions. 
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The DSE for the ghost propagator can be derived from the following operator identitj 
for the ghost kernel G: 

Gt = k-'S'^' + k-\r'^/rA'G^: (9.6) 

which follows from definition ( 17.5^ and Eq. (I7.2ip . Taking the expectation value of Eq. 
( 19. 6p . we find with the decomposition (I7.22p of the connected ghost-gluon vertex 

dj'^b = S'^^ + (T'^/y^DA^iTir'^ . (9.7) 

After contraction with 5°'^ and using the definition ( 17.250 for the form factor 7^ of the 
ghost-gluon vertex, Eq. ( 19.71) turns into 

, . / trTO'^rj; \ ^ dl . „ d 



n n 



In Fig. [9l the ghost Dyson-Schwinger equation ( I9.8p is depicted. It is equivalent to the 
exact ghost DSE in 3 + 1 dimensions. This is due to the fact that Eq. (|9.8p follows from the 



Since in I? = 1 + 1 pure Coulomb gauge the gauge field is constant, the convolution integral of 
the loop breaks down into a simple product in momentum space, see Eq. (|9.5I1 . 
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Alternatively, one can introduce ghost fields and proceed similarly as for the gluon propagator. 
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Fig. 10. Gribov regions of the pure Coulomb gauge, sketched as spherical shells defined by the 
magnitude of |A| = ki, k2, k^, k^, ... . The location of the fundamental modular region (FMR) 
is indicated by a fat line. 



operator identity (I9.6P and not from the details of the wave functionalliU The ghost DSE 



flQ.Sp can be solved for 7„ which confirms the relation fl7.26p found in section [7] for Nc = 2. 

In the derivation of the DSEs (19 .Sp and (|9.8p the integrated configuration space was set 
to be the first Gribov region fli of the pure Coulomb gauge, given by |A| < ^ = fci (see 
section HI). We now point at an important property of the Dyson-Schwinger equations. If 
the configuration space, here being the first Gribov regions Qi, is replaced by the union 
of the first two Gribov regions, Qi U f22, the DSEs (19.50 and (19.81) do not change. This is 
most readily seen for the ghost DSE (19. Sp which follows from the operator identity (19. 6p 
and therefore is not affected by the choice of the configuration space. The gluon DSE 
(19. 5p is derived from the path integral identity (19. ip which makes use of the fact that the 
Faddeev-Popov determinant jTp vanishes at the first Gribov horizon dfli. Changing the 
configuration space to Qi U Q2, the path integral identity (19. ip still holds true since by 
definition, the Gribov horizon is where the Faddeev-Popov determinant vanishes, 

Jp[A e d^n] =0 Vn . (9.9) 

We are therefore led in Q1UQ2 to the same gluon DSE as in Eq. (19. 5p . More generally speak- 
ing: Regardless of the configuration space, so long as it is a union of Gribov regions, Un ^n, 
the Dyson-Schwinger equations stay form invariant. This form invariance also applies to 
the DSEs in 3 + 1 dimensions. 

In Fig. [TOl the various Gribov regions fi„ of the pure Coulomb gauge are sketched. For 



On the other hand, the gluon DSE as it stands in (19.41) is only true for the actual constant wave 
functional. 



45 



definiteness, we put to our disposal the union 



Tng := ^^1 U U . ..^Na 



(9.10) 



of Ng Gribov regions. By the choice of Nq the configuration space is restricted to |A| < 
= k^G- Despite the form invariance of the DSEs with respect to T^q, their solution 
cannot be expected to be the same for all Nq- As a matter of fact, the gluon and ghost 
propagators strongly depend on Nq, as shown in the upcoming section. Thus, the set of 
DSEs does not have a unique solution, but for every Ng there exists (at least) one separate 
solution. The information on the configuration space (19.101) is missing in the set of DSEs 
and must be provided by subsidiary conditions, i.e. put in by hand. In this sense, the 
DSEs alone do not provide the full non-abelian quantum gauge theory. In the D = 3 + 1 
dimensional case, some approximation is used to solve the DSEs. Having obtained an 
approximative solution, there is no means of deciding on the value of Nq, i.e. whether this 
solution approximates the exact solution in the first Gribov region Fi = or rather in 
or any other T^^, remains unknown. The effect of truncating the set of DSEs in D = 1 + 1 
will be studied in section dH 



10 Many Gribov copies 

In this section, the configuration space is extended from the first Gribov region to a union 
^Ng (19.1 op of several Gribov regions, thus including many Gribov copies. Using the exact 
constant wave functional, the propagators, vertices and the colour Coulomb potential are 
calculated in T^^. Moreover, we extend the configuration space to Nq = oo, damping 
large Nq contributions with a Gaussian wave functional. This will illustrate the effect of 
insufficient gauge fixing on the infrared features of the theory. 



10.1 Exact Green functions 

The calculation of the Green functions in the extended configuration space T]y^, defined 
in Eq. (19.101) . is identical to the one in section [7] with the exception that the constant wave 
functional ( 16.2ip needs to be normalised differently. This can be accounted for by letting 




(10.1) 



in the expectation values (16.22p . 
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Fig. 11. Left: Exact ghost form factor, including Nq Gribov regions. Right: The form factor 7^ of 
the ghost-gluon vertex for various choices of Nq. 



For the gluon propagator in the pure Coulomb gauge we thus find from Eq. ( I7.2p 



(10.2) 



and it reduces to the result (17.41) for Nq = 1. Note that there is a strong dependence on 
the parameter Nq, i.e. a strong Gribov copy effect. 

The result for the ghost form factor dn in the pure Coulomb gauge follows from making the 
replacement (110. ip in the expectation value ( 17. lip . 



G) 



1 + 



4 1 



ttNo 



d-d 



■{p' sin^ ■§ 



:io.3) 



Obviously, the ghost form factor still approaches tree-level for the ultraviolet limit n — > 00. 
There are, however, substantial changes in the infrared. The result (110.31) for the exact 
ghost form factor dn can be seen in Fig. [TTJ In the first Gribov region, Nq = 1, the ghost 
form factor shows an infrared enhancement, as already shown in section [7l Including further 
gauge copies makes dn peak for intermediate momenta. This peak resembling a resonance 
appears at n = Nq where the momentum kn equals the radius of the configurations space 
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For n < Ng, the ghost form factor drops below tree-level, due to negative eigenvalues of 
the Faddeev-Popov operator. In the limit Nq 00, the lowest momentum mode dn=i 
approaches a definite value, 

hm d^iNG) = \. (10.4) 



It is interesting to note that the deviation of the ghost form factor dn from tree-level, when 



In the thermodynamic hmit, L — > 00, the peak of the ghost form factor in the first Gribov region 
is asymptotically at A; = 0. This corresponds to the "horizon condition" in Z? = 3 + 1 dimensions. 
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Fig. 12. Coulomb potential Vc{r) 



summed over all modes, gives a constant independent of Nq, 



3nN, 



■kNq 



■i?^ sin^ 



4 1 



ttNg 



G 



dd sin d{sm d — d cos d) 



(10.5) 

This is illustrated in Fig. [TT] and can be understood as being a consequence of the gluon 
propagator DSE. Solving the ghost DSE (19.81) for the ghost-gluon vertex 7„, 



7n 



kl {dn - 1) 



;i0.6) 



and plugging it into the gluon DSE (19.51) directly leads to Eq. (110.51) . 

Relation (110.61) holds for any SU{Nc) gauge group and, more importantly, for any choice 
Fatq of the configuration space. Having calculated the gluon propagator Da and the ghost 
form factor dn as functions of Nq, Eq. ( 110.61) gives the solution for the form factor 7„ of the 
ghost-gluon vertex. The result is shown in Fig. dH While within the first Gribov region, 
Ng = 1, there is hardly any deviation from tree-level, this deviation is quite pronounced 
for Ng > 1. Let us note here that an approximation of the proper ghost-gluon vertex 
by the tree-level vertex is good if and only if the configuration space is restricted to the 
first Gribov region Qi. Working with this approximation in solving the Dyson-Schwinger 
equations has an important effect on the propagators. This will be discussed in the next 
section. 



The form factor 0„ for the Coulomb propagator with the result ( 17.17p for Ng = 1 is found 
for general Ng by making the replacement (110.11) in Eq. (17.171) . By Fourier transformation, 
the Coulomb potential VG{NG,r) for external charges separated by r can be obtained in 
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the same manner as presented for Nq = 1 in section [3 The result is 



L-2r sin (2ttNgt'^ 

^Ng (10.8) 

and can be seen in Fig. [I2l While the string tension remains the same for all values of Nq 
(cf. Eqs. (18. 9p and p0.8l) ). visible effects occur for large ratios j. For Nq > 1, there are 
locally stable minima of the potential near r = -j, an unphysical gauge copy effect. 




10.2 Insufficient gauge fixing 



If lattice calculations use a gauge fixing, the common technique is to minimise a suitable 
functional along the gauge orbit. The restriction to the fundamental modular region is 
achieved only at the absolute minimum of this functional. If this absolute minimum cannot 
be reached exactly, extra gauge copies will alter the result. In order to mimic the situation 
that the infiuence of gauge copies cannot be strictly excluded but only suppressed, we 
choose here a Gaussian damping of the gauge field configurations, 



Af^exp(^-^Y Ar2 = ^(l + coth(a2)) , (10.9) 



extending the configuration space, see Eq. (19.10p . to Foo. Expectation values in the state 
(110. 9p are most readily obtained by replacing 



oo 



... ^ dde'^ ... (10.10) 



in the corresponding integrals. By adjustment of the (free) parameter a G R, it can be 
controlled how many Gribov copies have a considerable weight in an expectation value. A 
large value a will take along many Gribov copies while a small value localises the weight 
of field configurations around A = 0, suppressing Gribov copies. 

Let us calculate the constant gluon propagator Da f l7.ll) in the Gaussian wave functional 
(110. 9p as a function of a. Using the replacement pO.lOp in the expectation value (17. 2p gives 

DA{a) = Im^ (1)7^^ sin^^^^e-^/-^ = ^ ^ + ^^^) ' ^'"^-''^ 



Obviously, the larger the width a of the Gaussian (110. 9p . the larger D^(a). This agrees with 
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Fig. 13. Ghost form factor c?„(a) in the Gaussian wave functional (110. 9p with those choices 
ai, 02, 05, aio that yield the gluon propagator in the exact vacuum state with Nq = 1, 2, 5, 10. 
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a n 

Fig. 14. Left: The infrared enhancement, ^1(0;), is smaller than the physical value, di 1.37 for 
all values of a. Right: The ghost form factor in the state (jlO.Qh for the value a 1.40 that gives 
the strongest infrared enhancement, compared to the exact physical result. 
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the result (I10.2P of the exact wave functional where Da(Ng) rises quadratically with the 
number Nq of Gribov copies. One can now adjust a in Eq. (llO.lip such that it equals the 
result f ll0.2p for the exact vacuum state, given a specific Nq. For the values Nq = 1, 2, 5, 10, 
respectively, the corresponding values for a are 

ai^l.63, as -3.56, ~ 9.04 , aio^lS.l. (10.12) 

This procedure simulates the inclusion of Nq Gribov regions in the expectation values by 
a choice of the Gaussian wave functional (110.9^ . 

We now proceed to calculate the ghost form factor dn{a) in the state (llO.Op . The replace- 
ment pO.lOp in the integral (17. lip yields 

d^{a) = l + ^Af' L^sin^^ . .f .. e-'^'/^"'. (10.13) 
^ ^ 

The limits of a — and a ^ oo produce results that can be anticipated. As a 0, 
the Gaussian picks out the point A = from configuration space. Hence the tree-level 
behaviour of the ghost form factor appears: 

limc/„(a) = l, Vn. (10.14) 

a—^0 

The other extreme, a ^ oo, takes infinitely many Gribov copies into account and therefore 
must resemble the case Nq —>■ oo for the exact vacuum state. Indeed, for di{a), we find 

lim di(a) = - , (10.15) 

«^oo 3 

in agreement with Eq. (110. 4p . 

In Fig. [m the result (110.130 for (i„(a) is shown for the four values of a in Eq. (110. 12p 
which yield the exact gluon propagator for Nq = 1,2,5, 10. This should be compared to 
the ghost propagator in the exact wave functional for Nq = 1,2,5,10 in Fig. [TTl The 
effect visible in the exact wave functional, that taking more Gribov copies into account 
damps the infrared enhancement of dn and produces a spurious peak (here weakened) at 
intermediate momenta, can indeed be mimicked by the wave functional (110.90 with the 
appropriate Gaussian damping a. 

However, if one tries to quantitatively achieve the infrared enhancement of the physical 
solution, i.e. the exact dn within the first Gribov region fli, the wave functional ( 110. 9p 
fails. There exists no value for a such that the infrared enhancement of the exact physical 
solution, dn in Fig. [3l is realized. The value dn=i = 1.37 of the exact form factor ( 17. lip 
is larger than any choice of a can produce for the mode dn=i{cx) in the Gaussian wave 
functional (110. 9p . In Fig. [HI it is shown how di{a) varies with a. For the value a^ax where 
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Fig. 15. Coulomb potential Vc{a,r) in the Gaussian wave functional (|10.9I) with the "best" value 
(|10.16p for the width a. Also shown are the results in the exact ground state for the pure Coulomb 
gauge and the diagonal Coulomb gauge from Fig. [7l 



the infrared enhancement of the ghost propagator in the Gaussian wave functional (llO.Qp 
is the strongest, 

amax ~ 1.40 , d^iamax) ~ 1-20 , (10.16) 

it is seen how it still underestimates the physical result, di{amax) < 1-37. This indicates that 
if a lattice calculation is unable to exclude all Gribov copies, the genuine infrared physics 
cannot be described. In higher dimensions, this would mean that the infrared enhancement 
of the ghost form factor on the lattice is weaker than expected from continuum studies, an 
effect that is indeed observed [T9ll46l. 



Expecting that with the choice a = a^ax, most (though not all) of the infrared features can 
be carried alongmi we go on to compute other expectation values of interest. For instance, 
the Coulomb potential Vc{Q.,r) between two external static colour charges separated by 
r can be computed by taking the pure Coulomb gauge expectation value of the operator 
(18. 7p in the Gaussian wave functional (110. 9p . 
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This choice is reminiscent of the "horizon condition" in higher dimensions. 
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Fig. 16. Ghost-gluon vertex form factor 7,1 in the Gaussian wave functional (|10.9I) with the "best" 
value (110. 16|) for the width a. Also shown is the physical solution, see Fig. [U and the tree-level 
value. 

In the thermodynamic limit L — 00, the potential (110.171) approaches the same behaviour 
as Vc{r) in the exact ground state. It was already shown above that the Coulomb string 
tension ac, defined in the thermodynamic limit, turns out independent of the wave func- 
tional and the configuration space T^^. For the choice pO.161) of a, the potential Vc{r) in 
Eq. (110.171) is shown in Fig. [151 Varying a, the potential is seen to vary between the exact 
solutions of the pure Coulomb gauge and the diagonal Coulomb gauge. As a matter of fact, 
the result (I8.12P of the diagonal Coulomb gauge is reached in the limit a — *• 0, 

lhn\/c(«,r) = y - ^) = yc{r)U,on.i„ ■ (10.18) 

Recall that the limit a ^ turns the Gaussian wave functional (110.91) into a delta distri- 
bution, peaked at A = 0. The limit ( 110.181) therefore emphasises that in 1 + 1 dimensions, 
the Coulomb potential Vc{r) is a quantity that is independent of quantum fluctuations, 
when properly evaluated in the fundamental modular region. 

In the opposite limit, a — >■ cxo, the potential Vc{a,r) coincides with Vc{NG,r) in the exact 
vacuum state, see Eq. (110.81) . having taken the limit Nq 00. This is intuitively clear 
since the Gaussian state l\10M becomes constant for a — 00 (and not normalisable). 
Note, however, that the limits a — 00 and L ^ 00 are not interchangeable. In order to 
find the right string tension, one is to first evaluate the energy in a normalisable state 
{a < 00), take the thermodynamic limit L — ^ 00, and then determine the string tension as 
a function of a. In the opposite order of limits, one finds a different result for ac- 

Finally, we calculate the ghost-gluon vertex in the state pO.91) . It can be found using the 
expression (110. 6p which holds for any wave functional (and any T^^); in the evaluation we 
simply use the results (llO.lip and (110.131) for the gluon and ghost propagators. The result 
is shown in Fig. [161 The deviation of 7„ from the exact physical result is small, although it 
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is enhanced for the infrared modes. Anyhow, the tree-level vertex is a better approximation 
of the exact ghost-gluon vertex than the one in the Gaussian wave functional (110.91) . 



11 Truncation effects 



In higher dimensions, it is not possible to obtain exact nonperturbative expressions for the 
Green functions. A common approximation that enables us to find the infrared asymptotic 
solutions of the Dyson-Schwinger equations is to render the ghost-gluon vertex tree-level. 
We will make this approximation here in 1 + 1 dimensional pure Coulomb gauge and inves- 
tigate the effect on the propagators. To be explicit, we set 

r;: - , ^ 7n = i • (ii.i) 

Recall it was shown in the preceding sections that despite the form invariance of the exact 
DSEs with respect to the configuration space ^Ng^ the true Green functions (calculated 
with the exact vacuum wave functional) do depend on Ftv^. With the approximation (111.11) . 
the DSEs for the gluon and ghost propagators, Eqs. (19. Sp and (19.81) . turn into 

d^ = l + N,DA^, Df = N,T.%- (11-2) 

Diagrammatically, these equations are depicted in Figs. [8] and [H with the blobs replaced 
by dots. We note that due to the approximation ( 111.11) . this set of equations is closed 
and can be solved. The crucial point is to realize that the solution will no longer depend 
on the choice of the configuration space Tj^^. Since the ghost-gluon vertex was chosen 
to be independent oiV^^, see Eq. (Ill.ip , the gluon and ghost propagators are now also 
independent of FArg. The original exact set of DSEs holds within any of the Gribov regions. 
However, solving the approximated set of DSEs, the information is lost in which Gribov 
region the Green functions are evaluated in. This problem also occurs in the infrared ghost 
dominance model of the D = 3 + 1 theory. 

Let us now solve Eq. pi.2p explicitly. Since it is quadratic in (i„, we find two solutions. 



1 ± /I _ 4iY^ 

dn = . (11.3) 

In order for the limit (i„ ^ 1 for A;„ = ^ cxd to be fulfilled, we need the lower sign. 
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Fig. 17. Left: Comparison of the exact gluon propagator to the one from the DSEs with a tree-level 
ghost-gluon vertex, shown as a function of Nq. Right: Comparison of the exact ghost form factor 
dn in the first Gribov region [Nq = 1) to one calculated from the truncated DSEs, both shown 
as functions of the momentum mode n. 

Plugging this solution into the second DSE in Eq. 011.2p . we find 

2 



1 



\ 



1 



[Tm] 



DaN^L'^ ^ 7.74 



:il.4) 



The numerical solution yields with = 2 for the gluon propagator 

1 



3.87 



L2 



(11.5) 



This value can be plugged into Eq. (Ill.Sp . with the lower sign, to immediately get the 
result for the ghost form factor. 



The gluon propagator ( 111.5p is evidently independent of the Gribov region. One can choose 
with any Nq, the approximation pi. II) will always yield the gluon propagator given 
by Eq. (111.5p . On the other hand, we know the exact result for the gluon propagator with a 
given value of Nq, see Eq. (110. 2p . The latter can now be compared with the approximative 
result (lll.5p . for each No- 

In Fig. [Ill it is seen how for large Nq, the exact and approximative results differ dramat- 
ically. However, for Nq = 1, i.e. in the first Gribov region, there is a good agreement. The 
tree-level approximation for the ghost-gluon vertex yields a value for the gluon propagator 
that is very close to the exact result in the first Gribov region. The same occurs for the 
ghost propagator. As can be seen in on the right panel of Fig. [ITl the approximative (with 
7n = 1) DSE result for the ghost form factor dn hardly deviates from the exact result for 
Nc = 1. At the same time, the approximative and exact solutions of dn for other values of 
Nc disagree both quantitatively and qualitatively, cf. Fig. [TTl Also shown in Fig. [TT] is the 
ghost-gluon vertex in the exact vacuum state for different values of Nq. One can realize 
that the tree-level approximation (111.10 is only for Nq = 1 a good one. 
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We infer that while the exact Dyson-Schwinger equations are form invariant with respect to 
the configuration space T^^, the approximation of the ghost-gluon vertex by its tree-level 
value effectively puts the approximative solution of the propagators into the first Gribov 
region Fi = This supports the approach of solving DSEs in D = 3 + 1, truncated by 
means of a tree-level ghost-gluon vertex. 



12 Variational approach 



In order to test the variational approach to Yang-Mills theory in 3 + 1 dimensions, the 
same ansatz for the wave functional as in Ref. [8] is here applied to solve the Yang-Mills 
Schrodinger equation in D = 1 + 1 in the pure Coulomb gauge. The variational calculation 
is performed using the same approximations as in Ref. [8j. 

The variational wave functional is given by 

^(A) = , ^ ^(A) , ^{A) = ATe^s^""^" := {A\uj) , (12.1) 
Jp{A) 



where u is a variational parameter determined by minimising the vacuum energy density. 
The ansatz (112.11) is mainly motivated by simplicity. It removes the Faddeev-Popov deter- 
minant J^p from the integration measure in the expectation values of the pure Coulomb 
gauge, 

(*|C)(A)|^) = J VAJp{A)0{A)\^{A)\^ = J VAO{A)\^{A)\^ =: {0{A))^ . (12.2) 

Too Too 

and thus allows for an immediate application of Wick's theorem. Here and in the following, 
{. . .)u) denotes the expectation value in the state \u) (|12.1I) . for which the scalar product 
is defined with the flat integration measure VA = HdA"^. 

a 

The conflguration space in the expectation value 012. 2p is not restricted to the first Gribov 
region fli but is extended to the union Fqo of all Gribov regions, see Eq. (IQ.lOp . This is 
motivated from the -D = 3 + 1 case where only little is known about the Gribov horizon 



and a restriction to fli is technically cumbersome[i£J In section 110.21 we have seen that 
within Foo, any Gaussian damping of gauge copies will fail to recover the exact infrared 
behaviour of the ghost form factor. Note, however, that the wave functional in Eq. p2.1l) is 
supplemented by a Faddeev-Popov determinant. We will show below that with appropriate 
approximations, the correct infrared behaviour of the ghost form factor can thus still be 
maintained. 



The so-called Gribov-Zwanziger action is used in a few studies [l^ to realize a restriction to 

ni. 
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The normalisation constant M in the ansatz (112.10 for the wave functional is chosen such 
that (1)^ = 1 and is given by 



In AT 



,7r 



Thus, the static (equal-time) gluon propagator reads 



(12.3) 



(12.4) 



In D = 3 + 1 the kernel u has the meaning of the gluon energy. In the present 1 + 1 
dimensional case u has dimension mass~'^ and is required by normalisability of the wave 
functional to be positive, > 0. Let us emphasise that whatever the variational principle 
yields for cu, it will determine the gluon propagator Da by Eq. p2.4p . 

The vacuum energy E{uj) is calculated by taking the expectation value of the Yang-Mills 
Hamiltonian in the pure Coulomb gauge f l5.20p in the absence of external charges. After a 
partial integration, Eiuj) yields 

E{uj) = (^|//|^) = J VAJp{A)^*{A)H'^{A) 



VA (ni*)* (ni^) 



where 



Using 



we find 



m = jp/'nij^'/' = ni - ^ (ni in jp) 



1 

2 



ni^ = - 

1j 



■ 1 

2 I dA'^ 



(12.5) 

(12.6) 
(12.7) 



9_ 
2L 

2L 



2 rfv4« 



^2 (^A^A"")^ + uj{A 



dA" 



1 I ( dlYiJp 
^ 4 \ I 



;i2.8) 



Following Ref. [8] we will explicitly calculate the first two terms and then find the last term 
by completing the result to a total square, which is correct up to two loops. The first term 
in Eq. (I12.8P can obviously be expressed by the gluon propagator (112. 4p . For the second 
term, we use the abbreviation 



X 



ab 



d\njp 

dA^ 



(12.9) 
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which was referred to as the "curvature" in Ref. [8]. Using the definition (17.221) of the proper 
ghost-gluon vertex and its form factor 7„ (17.250 . the curvature can be written aa ^^1 

= + c {A^ Tr G r°/)^ = ^ E 7n § tr (r^r^'' 



= VE^"^^'^'=^^^'^' (12.10) 

and its diagonal elements define the scalar curvature x- Performing the quadratic comple- 
tion, the expression p2.8p for the vacuum energy can be cast into the form 

E{oo) = g' ^j—^- ^ (12.11) 

and it is obviously minimised for the choice 

uj = x (12.12) 

of the variational kernel u. Equation (|12.12p is called the gap equation and it gives rise to 
an infrared divergent gluon energy uj{k) in 3 + 1 dimensions [8|27| . 

The gap equation ( 112.121) states that the gluon propagator ( 112. 4p can be related to the 
curvature x and we may use the definition p2.9l) of x°'^ to calculate the gluon propagator 
exactly. However, approximations have been made and — more importantly — the configu- 
ration space was not properly restricted to Fi. Let us look at the task of determining the 
solution for the gluon and ghost propagators differently. We note that plugging the gap 
equation ( fT2T2D into Eq. ( fl^ yields with Eq. STTm 

Dt = NcY.yn^- (12.13) 

Turning back to Eq. ( 19.51) . we recognise that relation p2.13p is identical to the gluon 
propagator DSE in the exact vacuum state. Moreover, we can use the ghost propagator DSE 
(19. 8p from the exact vacuum state since it follows from an operator identity, independent of 
the wave functional or configuration space. The set of equations the variational calculation 
above resulted in is equivalent to the set of Dyson-Schwinger equations derived in the 
exact vacuum state. The difference is that here the expectation values (i.e. Da, 7n) 
are evaluated in the state ( 112. ip and in the configuration space Too, yielding different 
results. It was shown in section [10] that choosing a set of several {Nq > 1) Gribov regions 
for Ftvc results in drastic changes for the Green functions. However, if we use the tree- 
level approximation for the ghost-gluon vertex, it is clear from the discussion in section 



The definition (|12.9p of the curvature is equivalent to the one in Ref. [8] within the wave 
functional p2.1|) . The same holds for the proper ghost-gluon vertex. 
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[TT] that the solution to the DSEs so obtained is very close to the exact solution. With 
a conspiracy of approximations, namely the quadratic completion in Eq. (112. Sp and the 
vertex approximation 7„ = 1, the variational state (112.11) in Too yields the same propagators 
as the exact vacuum state (16.211) in the first Gribov region Fi =Qi. 

Identifying the variational wave functional 012. ip for the solution u = x with the exact 
wave functional = const implies that the Gaussian must cancel the Faddeev-Popov 
determinant jTp, 

Jp exp • (12.14) 

In Ref. [I^ it was shown that up to two-loop order in the energy the replacement p2.14p 
is exact and thus results in the correct DSEs. 

_ 1 

In D = 3 + 1 both u and x are momentum dependent and the cancellation of J'p ^ against 
the Gaussian in the wave functional is obtained in the infrared limit k ^ only. We thus 
observe that in the infrared limit the wave functional in D = 3 + 1 reduces to the exact 
wave functional in D = 1 + 1 dimensions. As discussed in Ref. [10] the constant wave 
functional does not constrain the infrared modes of the gauge field and thus describes a 
stochastic vacuum where the infrared modes can arbitrarily fiuctuate. 

As shown in Ref. ^U\, the cancellation of Gaussian and Faddeev-Popov determinant per- 
sists in D = 3 + 1 in the infrared even if the more general ansatz is used, 

^(A) = J7"(A)§(A) . (12.15) 

In this state with a arbitrary exponent a of the Faddeev-Popov determinant, the gluon 
propagator becomes 

{A'^A'') = 5''\2^)-^ , (12.16) 

where 

u = uj-{2a-l)x- (12.17) 
The gap equation is the same as above, see Eq. (112.121) . except that u is replaced by uj, 

u = X. (12.18) 

We therefore find from 012. 17p 

uj = 2ax. (12.19) 

For a = I we recover, of course, the previous result (112. 12p . while a = yields u = and 
the variational wave functional (112. 15p becomes the exact one 

"^{A) =Af = canst. (12.20) 

and thus yields also the exact results for the propagators, provided the range of the field 
A is properly restricted to the first Gribov region. 
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Finally, let us turn to the Coulomb form factor /„ which measures the deviation of the 
Coulomb propagator {G{—d'^)G) from the factorised form {G) (—5^) {G) and was calculated 
in the exact wave functional at the end of section [731 In 3 + 1 dimensions, the (momentum- 
dependent) form factor f{k) is set to unity since it fails to satisfy the corresponding integral 
equation within the approximations made |48| . The form factor f{k) requires a higher-order 
calculation (as pointed out in Ref. P|3T] ). and for this reason it is investigated here in 1 + 1 
dimensions where approximations are not necessary. 

The integral equation that is derived for the Coulomb form factor fn follows from the 
identity [49] 

FiA) = GiA)i-d')G{A) = ^(gGigA)) . (12.21) 
Here, we have scaled the gauge field by the coupling constant g, 

A = gA (12.22) 

so that with 

G-\gA) = ~d^-gAd (12.23) 
we can derive Eq. (112.211) by differentiation. Following Ref. |2j, in the variational approach 
[8j the vacuum expectation value of the relation (|12.2ip was taken, thereby ignoring the 
implicit ^f-dependence of the wave functional to obtain the approximative relation 

{F)^^^ig{G)) (12.24) 

which is the so-called Swift relation p]. In the present 1 + 1 dimensional case the exact 
vacuum wave functional is independent of g but the Faddeev-Popov determinant J'p{gA) 
in the integration measure is ^f-dependent and this ^f-dependence is ignored in Eq. f ll2.24p . 
Expressing {G) and (F) in terms of the gho st and Coulomb form factors, dn (17.81) and /„ 
(EIHj), the Swift relation ( 112.241 ) become^" 
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fn ~ d~' — {gd^) = d-'-g—d-' . (12.25) 



Using the inverse form of the DSE (19.81) for the ghost form factor, 



djn 



d-^ = l-N,^nDA^ , (12.26) 



an integral equation for /„ can be found, 



d^ f 



/„^l + iV,7„D^^. (12.27) 



^® This relation differs from the one given in Ref. [8], where an extra factor of g was included in 
the ghost form factor. 
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Fig. 18. Coulomb form factor /„, defined by Eq. (IT.lSp . in the exact vacuum state depending on 
the number Nq of included Gribov regions. 

Note that this integral equation is not exact, due to the Swift approximation (I12.24p . In the 
exact wave functional of the 1 + 1 dimensional pure Coulomb gauge, the ghost propagator 
is independent of g ^ 
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d J^^VAG{gA)DetG~\gA) 
dg J^^VADetG-^igA) 
d J^^VAGiA)DetG-\A) 
dg !^^VADetG-\A) 







;i2.28) 



and hence the ghost form factor also is, ddn/dg = 0. The Swift relation (112.251) thus 
simplifies to 

fn = d-' . (12.29) 
This relation implies that /„ is infrared suppressed if dn is infrared enhanced, which is in 
contradiction to the true behaviour of fn and dn, see Figs. [3] and El Even more directly, the 
contradiction can be seen by plugging Eq. (112.291) into the approximative integral equation 
(112.271) and comparing to the exact ghost form factor DSE (112. 26p . 

The contradiction arises from the approximation made to arrive at the Swift relation 
(112.241) . In the infrared limit of the D = 3 + 1 theory which is correctly described by 
the D = 1 + 1 wave functional [lO], the Swift relation therefore must lead to inconsisten- 
cies [3T1I481. 



Recent lattice calculations in = 3 + 1 of the Coulomb form factor f{k) seem to indicate 
that in the infrared f{k) is enhanced. This enhancement gets weaker the better the Coulomb 
gauge is fixed. In 1 + 1 dimensions, the form factor /„ can be calculated exactly, see section 
I7.5[ The effect of including several Gribov regions, i.e. choosing rArQ>i as done in section 



The same is true for the ghost propagator in 3 + 1 dimensions if we consider the stochastic 
vacuum "^[A] = const of the Coulomb gauge Hamiltonian approach or the ghost dominance 
model Sym = in the Landau gauge to study the infrared asymptotics. 
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110.11 leads to an interesting observation. In Fig. [THl it is shown how the exact result for /„ 
varies with the number No of included Gribov regions. The larger Nq, the more pronounced 
a (spurious) infrared enhancement. This gauge copy effect is in agreement with the findings 
on the lattice in 3 + 1 dimensions and indicates that in order to get the exact result for 
/(/c), gauge fixing on the lattice has to be performed very carefully. 



13 Summary and Conclusions 

In this paper, we have considered 1 + 1 dimensional SU (2) Yang-Mills theory in canonical 
quantisation in the pure Coulomb gauge as a testing ground for Yang-Mills theory studies 
in higher dimensions. The investigations were carried out in the pure Coulomb gauge and in 
the diagonal Coulomb gauge, where the residual global gauge invariance, left unfixed in the 
pure Coulomb gauge, is fixed by diagonalising the constant spatial gauge field. Although the 
two gauges differ only by a global SU{2)/U{1) gauge-fixing constraint, they have different 
Faddeev-Popov determinants due to additional zero modes of the Faddeev-Popov kernel in 
the pure Coulomb gauge. While the pure Coulomb gauge is perfectly suitable for perturba- 
tion theory, the diagonal Coulomb gauge is ill-defined for the perturbative vacuum A = 0, 
for which the Faddeev-Popov determinant vanishes. The occurrence of such gauge-fixing 
defects is a characteristic feature of so-called abelian gauges the diagonal Coulomb gauge 
belongs to. In higher dimensions gauge-fixing defects of abelian gauges manifest themselves 
as magnetic monopoles in the corresponding abelian projection [50j (see also Ref. [51i|). 

We have explicitly demonstrated that the Faddeev-Popov method does not require com- 
plete gauge fixing but works for any partial gauge fixing, provided that the zero modes of 
the Faddeev-Popov kernel arising from the residual gauge symmetry (left unfixed by the 
partial gauge fixing) are properly treated. In the resolution of Gauss' law, these zero modes 
give rise to residual constraints on the wave functional, which express the invariance of the 
wave functional under the residual gauge symmetry: The Noether charges corresponding 
to these residual symmetries must vanish in physical states. The constraints on the wave 
functionals arising from the residual unfixed gauge symmetry exist also in higher dimen- 
sions but have not been explicitly identified so far, except for space-independent gauge 
transformations [35j. They also naturally emerge in the functional integral approach in the 
so-called first order formalism where the temporal gauge field can be explicitly integrated 
to leave a 5- functional, which enforces Gauss' law [T4j. In the pure Coulomb gauge, the 
Gauss' law constraint can be worked out analogously to the Hamiltonian approach and 
the 5-functional can be used to integrate out the longitudinal components of the momen- 
tum field. When the resolution of Gauss' law is properly done, i.e. the zero modes of the 
Faddeev-Popov kernel properly treated, from the ^-functional some ordinary (5-function 
survives, which precisely enforces the vanishing of the Noether charges corresponding to 
the residual unfixed gauge symmetries [35] • 
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The exact spectrum of the Yang-Mills Hamiltonian was obtained within both the diagonal 
and the pure Coulomb gauge, having implemented the constraints on the wave functional 
arising in the resolution of Gauss' law from the zero modes of the Faddeev-Popov kernel. 
In the thermodynamic limit, we recovered the well-known spectrum that leaves only the 
vacuum state at zero energy, freezing out all excited states. The exact vacuum state was 
used to calculate the ghost and gluon propagators, the ghost-gluon vertex and the static 
colour Coulomb potential. We compared the results in the pure Coulomb gauge restricted to 
the first Gribov region to those in the diagonal Coulomb gauge restricted to the fundamental 
modular region. For the propagators, the colour trace was found to be left invariant when 
transforming from the pure Coulomb gauge to the diagonal Coulomb gauge. We found that 
the ghost propagator is infrared enhanced, in agreement with the horizon condition widely 
used in Dyson-Schwinger studies of D = 3 + 1. This infrared enhancement is the strongest 
when the configuration space is properly restricted to the first Gribov region. We studied 
the effect of including several Gribov copies, either by extending the configuration space 
to a union of Gribov regions, or by using all Gribov regions with a Gaussian damping. 
It was seen that the quantitative infrared enhancement of the ghost propagator cannot 
be realized by any of these calculations that include gauge copies from outside the first 
Gribov region. This indicates that lattice calculations of the ghost propagator require a 
very accurate gauge fixing and explains the shortcomings of the infrared enhancement of 
the ghost propagator on the lattice when compared to continuum studies [T9f46] . 

The Coulomb string tension yielded the same results for both gauges which is a fortunate 
result, considering that in D = 3 + 1 calculations the gauge is not completely fixed. The 
quantitative result of the static colour Coulomb potential, away from the thermodynamic 
limit, differs for both gauges. The result of the pure Coulomb gauge can be sort of artificially 
deformed into the result in the diagonal Coulomb gauge by suppressing Gribov copies with 
a Gaussian wave functional of width zero, as discussed in section [TOl The investigations 
showed that on x R, the Coulomb string tension arises from the abelian part of the 
Coulomb interaction and actually is identical to the string tension of the abelian theory, 
thus providing an upper bound of the gauge invariant string tension [44]. The effect of 
gauge copies on the static colour Coulomb potential was studied by taking several Gribov 
regions into account, and resulted in spurious locally stable minima for large separations 
of external colour charges. 

The Dyson-Schwinger equations for the propagators and vertices in the pure Coulomb 
gauge were derived. It was shown that the exact solution within the first Gribov region 
satisfies the Dyson-Schwinger equations, but that these solutions are not the only ones. 
Changing the configuration space from any union of Gribov regions to another leaves the 
Dyson-Schwinger equations form invariant. This persists in the D = 3+1 case. Therefore, it 
is legitimate to ask: In which union of Gribov regions are the Green functions given when 
the set of Dyson-Schwinger equations is solved by means of a truncation? We found in 
D = 1 + 1 that choosing the ghost-gluon vertex at tree-level, effectively puts the solutions 
for the propagators and vertices into the first Gribov region. As for the colour traces 
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of the propagators, even the result within the fundamental modular region is attained. 
The tree-level ghost-gluon vertex approximation for solving Dyson-Schwinger equations — 
which was advocated by several investigations before |8|52f53f54| — thus receives further 
strong support. 

The variational approach to Coulomb gauge Yang-Mills theory in = 3 + 1 dimensions 
[8] integrates over all Gribov regions for technical reasons. It was shown in section [12] by 
using the same ansatz for the vacuum wave functional in D = 1 + 1, that with the appro- 
priate approximation (quadratic completion in kinetic energy expression) the variational 
principle yields a set of Dyson-Schwinger equations that is the exact one. With the tree- 
level ghost-gluon vertex, the exact solution within the first Gribov region is thus very well 
approximated. It can be expected that the infrared limit of the D = 3 + 1 theory, which is 
described by a stochastic wave functional as in D = 1 + 1, is thus also well-approximated. 
Furthermore, we found that the Coulomb form factor that is necessary for the calculation 
of the static colour Coulomb potential is not too far from tree-level in the exact D = 1 + 1 
calculation and that the inclusion of many Gribov copies simulate a spurious infrared en- 
hancement of the Coulomb form factor. We infer there are no indications that the choice of 
a trivial Coulomb form factor is worse than any other approximation made in the D = 3 + 1 
calculations. 
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Appendix A SU{2) colour rotations 



The unitary matrix U, which rotates the colour vector n(^, 0) into the 3-direction (and 
thus, in particular, diagonalises the gauge field (I2.19p ) is defined by 

u^h{e,<p)-Tu = n. (A.l) 

This matrix is defined up to an abelian gauge transformation U Uuj , oj = explipT^) e 
U{1) C SU{2), i.e. it is defined on the coset SU{2)/U{1). The adjoint representation U, 
defined by 

U^TaU = U^^Tb (A.2) 
is related to the fundamental representation U by 

Uab = -2 tr (u^TaUn) , tr (T^T^) = -^6ab • (A.3) 

Equation (IA.3I) is most easily proved by Taylor expanding U in terms of 9 = QaTa and 
using 

[r,,e] = e"^r,. (a.4) 

Furthermore, since (1A.2I1 is based only on the algebra of the generators it is valid in any 
representation, in particular the adjoint representation (where W = U'^), 

U^tif = Uabfb . (A.5) 

The matrix U can be realized by 

[/(^, 0) = e^^^ ■ T , T = -V, (A.6) 

where 

60 = — sin ei + cos 62 (A. 7) 

is the unit vector in the direction of the azimuthal angle 0. The matrix U ( 1A.6I) can be 
alternatively expressed in terms of Euler angles as 

U{e, 0) = e^^^e^^^ . (A.8) 

Equations ( lA.ll) . (1A.6I 1. ( 1A.7I1 are valid in any representation of SU{2) and thus also in the 
adjoint representation 

U{e, 0) := e^""^ ■ * = e^^^e''^^ . (A.9) 

From the defining equations (lA.ll) . (IA.3I) and = also follows that the components 
of the unit colour vector n{9, 0) are given by 

n%e,(i)) = Ua3{e,<p). (A.io) 



65 



It is convenient to use the bracket notation 

IJab = {a\U\b) (A.ll) 
and to express JJ in the basis of the eigenvectors \a = 0, ±1) of the spin 1 operators 



(A.12) 



satisfying 



5» = 1(1 + 1)|(t) 
5*310-) = o\a) . 

The transition elements 

(a|a) =: e^ 

are the Cartesian components of the spherical unit vectors (in colour space) 



(A.13) 
(A.14) 



ea=i 



1 1'^ 



6^=0 



^0^ 



(A.16) 



Here, Greek letters a, r, . . . denote spherical colour components {1,0,-1}, while Latin 
letters a,b,... denote the Cartesian colour components {1,2,3}. The matrix elements of 
the adjoint representation U (|A.9l) in the spherical basis 



(A.16) 



(a|[/(e,0)|a') = {a\a){a\U\b){b\a') = ef f/„,(e, 0)e^ 
are related to the Wigner D-function by 

{a\u{e,<p)W') = Di:,\<p,e,o). 

Using ((t|3) = e^* = S^o, the colour unit vector ( lA.lOl ) can be expressed as 

n'^i9,<P) = {a\a){a\U\r){r\3) = e:DU<P,9,0) . (A.18) 



(A.17) 



Appendix B Explicit resolution of Gauss' law 



To identify nM(a;) we Fourier expand the periodic gauge field A{x + L) = A{x) 



1 27r'n 
A{x) = -Y.e''-^A{n), = , neZ 



(B.l) 
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The inverse transformation reads 

L 



A{n) = / dxe-'^-^Aix) (B.2) 







and the continuum Umit L — > oo is obtained by the replacement 



For later use we also quote the completeness and orthogonality relations 

1 1 r 

^(^) = L^ ' = L y e^C^'-^")- , (B.4) 

where 5{x) denotes the periodic (5-function, satisfying 

5{x + L)=5{x). (B.5) 
From f IB.ll) we find for the momentum operator 

n"(a;) = = T E e'^^"" .,,"^, , . (B.6) 
B.l Pure Coulomb gauge 

In momentum space the pure Coulomb gauge (12.16p reads 

A{n) = 5„,oA(0) , (B.7) 

where 

A(0) = -j dxA{x) =■ A (B.8) 



is the constant part of the gauge field, which is left after gauge fixing. From Eq. ( 1B.6I) 
we read off the transversal (x-independent) and longitudinal (x-dependent) parts of the 
momentum operator to be given by 

ni = , nt](x) = T E e^^^-'-r/rr • (B-9) 

In the pure Coulomb gauge, where the degrees of freedom are A'^^^''^'^, the charge of the 
gauge bosons ( 15.6p is space independent but non-zero. Using (cf. Eq. ( 13.26^ ) 

in^b^-ikux ^ e-ik^^-^{^a\lf\a)Xn,a{(r\U^\b) , (B.IO) 



67 



we obtain 

^b'^'Ii\{x) = W e-^'-^{a\U\a)KA^\U^\b)-T^, ■ (B.ll) 

Inserting this relation into Gauss' law f l5.4p and multiplying the resulting equation by e*^™^, 
and integrating over x thereby using Eq. (IB.4P we obtain 

d f 

^ 5^^MUW)KM\U^\h)-^^^^{A) = I j dxe^'-^pUx)^{A) . (B.12) 



For m = the l.h.s. of Eq. (IB.12I1 vanishes and we find that the wave functional has to 
satisfy the following constraint 

L 

Q'^^iA)^ JdxpU^)^{A) = 0. (B.13) 



For m 7^ the summation over n on the l.h.s. collapses to the term m = n and Eq. (IB.12P 
becomes after multiplying it by {a\U'^\a) and summing over a 



L 



\rn,MU^\^)-^^^^^{A) = z{a\U^\a) J dx e^''-^ p^Jx)^ (A) . (B.14) 
Multiplying Eq. flBlil) by le-'''"^y {c\U\a) X;^]^ and summing over m ^ and cr we get 

= ^7 E Ee-^''"' / dxe^'-^c\U\a)\;^\{a\U''\a)pUx)^>{A) (B.15) 

L 

= ijdxY. E^m,.(l/)A-;.<.(x)p»„,(x)^(A) , (B.16) 

where we have used the explicit form of the eigenfunctions ^-(x) (13.281) of the covariant 
derivative il)"-^ . Note if the modes m = 0,cr = ±1 were included in Eq. (IB. 161) . the sum 
would produce the inverse kernel 

{yMm-'\a.x) = j:'vlJy)Klv^A^) , (B.17) 

where the prime indicates that the mode m = cr = is excluded (while m = 0, o" = ±1 is 
included). 
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It is now straightforward to calculate the Coulomb Hamiltonian He defined by Eq. (15.91) . 
With Eq. (IB.16P we obtain after straightforward manipulations 

2 

Hc = yI dxdypl,{x)F^\x,y)pUy) , (B.18) 



where 

is the so-called Coulomb kernel. Let us stress the mode n = 0, a = ±1 is here not included 
although it is not a zero mode A„=o,cr=±i 7^ 0. Since this mode is also excluded from the 
ghost kernel 

G^'{x,y) = {x,a\G-'\y,b) = {x,a\{-Dd)-'\y,b) = EE<<.(^) (A;„A„,.)-^ 

(B.20) 

the Coulomb kernel ( IB. 191) can be represented as 

F'^\x,y) = {x,a\{-bdr\-d'){~Dd)-'\y,b) , (B.21) 

which is the usual representation. In the abelian case, this kernel reduces to the usual 
Coulomb potential. In the non-abelian theory, this is a dynamical object depending on the 
field variables via the covariant derivative. 



The above derivation of He has shown that Gauss' law in the pure Coulomb gauge does not 
only give rise to the Coulomb Hamiltonian He ( IB. 181) but in addition yields the constraint 
(IB. 131) on the wave functional. This constraint arises from the zero modes of the Faddeev- 
Popov kernel, which are a consequence of the incomplete gauge fixing. In a complete gauge 
fixing such residual constraints would not arise . 

Due to the fact that constant modes kn = are excluded from the Coulomb kernel, the 
dynamical charge of the gauge bosons, pg, being space-independent, drops out from the 

'^^ When the constraint (IB.13j) is obeyed by the wave functional the mode m = 0,a = ±1 can 
safely be included in the Coulomb kernel (IB.IOP since it does not contribute when the Coulomb 
Hamiltonian acts on the wave functional. Thus with the constraint (|B.13I) satisfied we can use the 
alternative kernel 

F-'ix,y) = Y,'{x,a\U\n,a)X~f,{n,a\U^\y,b) = {x,a\{iDr^\y,b) (B.22) 

n,cr 

in the Coulomb Hamiltonian. It is precisely this kernel (but with A restricted to the hyperplane 
of the diagonal Coulomb gauge) which arises as "Coulomb kernel" in the diagonal Coulomb gauge 
derived in the next subsection. 
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Coulomb term (IB.lSp . In fact, with the exphcit form of eigenfunctions (p'^^{x) (I3.28P we 
have 

/ dxF''\x, y) = \ fdxY, e'^-^''~y^ F^^ = ^ 6n,oF^^ = . (B.23) 
Thus we can replace in He ( 1B.18I ) the total charge Pj^t by the external charge p", 

2 

Hc = yI dxdyp\x)F'^\x, y)p\y) , (B.24) 

and in the absence of external charges p"(x) = the Yang-Mills Hamiltonian reduces to 
the transversal part (IG.lSp . 



B.2 Diagonal Coulomb gauge 



In the diagonal Coulomb gauge ( 12.20p the remaining physical degree of freedom of the 
gauge field is ^ = A^{n = 0) and the corresponding physical momentum reads 

ni = *-ni,ni = i-^. (b.2b) 

We will keep here the same notation as in the pure Coulomb gauge and denote the remaining 
unphysical part of the momentum operator by 

nj] = n° - n;^ . (b.26) 

This part is given here by 

nil = T E 'e~'''""-71^ ' (B-27) 

II L ^ idA'{n) ^ ' 

where the prime indicates that the term n = is excluded for the generator of the Cartan 
algebra a = Oq = 3 only. Note that contrary to the pure Coulomb gauge, the "transverse" 
components d/dA"-^°'{0) belonging to the generators a = a of the coset SU{N)/U{1)^~^ 
are here parts of IIn. 

With the explicit form of the gauge-fixed field fl2.20p and the corresponding momentum 
Il± (IB. 251) one notices that the charge of gauge boson (15.61) vanishes in this case, 



a ^ _^a^Yib^ ^ -if ni = -r^'^AlUl = . (B.28) 



Inserting the explicit form of nij, given by Eq. (IB. 2711 . into Gauss' law (15.41) and using 
furthermore 

■^a6j^3y^|g-ifc„x _ J2(^a\a)Xn,a{(T\b) , (B.29) 
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Gauss' law becomes 



^E'e"^'"^'E(«l^)V.(^l&)^^^*(^) =^P''W^(^) • (B.30) 



Multiplying this equation by e*'^"'^(cr|a), integrating over x and summing over a, we obtain 

L 

idA^{n) 



T.'<n^K,Mb)-^^^^'^{A) = i{a\a) I f/xe*'=-V(a;)*(^) • (B.31) 



Recall that the prime indicates that the term n = is excluded from the sum for 6 = 3. 
Since {a\b = 3) = e^* = S„o on the l.h.s. the term tt, = is excluded for a = 0. Thus for 
m = a = the l.h.s. vanishes and we find the following constraint on the wave functional 

L 

Q^^ = J dxp^{x)^ = , (B.32) 



which should be compared with the constraint f lB.131 1 in the pure Coulomb gauge. Since 
in the present case the charge of the gauge bosons vanishes, Eq. (IB. 321) is the restriction 
of the constraint (IB. 131) to the charge of the Cartan subgroup. For m = 0, cr 7^ and for 
m 7^ 0, (T-arbitrary, Eq. (IB. 3111 becomes 

d 

^-^A^\^) ,dA\mf ^^^ = ^{(y\o) I dxe'''-^p\x)^{A) . (B.33) 

Note that since (a 7^ 0|6 = 3) = the summation over 6 is for m = 0, a 7^ restricted to 
6=1,2. Multiplying the last equation by (c|cr)A~^^ and summing over a we obtain 



d 



idA'^im 



1^ 

-^{A) = / dxe^'^-^p'^{x)^{A) , (B.34) 



where the prime indicates again that the term a = is excluded for m = 0. Multiplying 
this equation by e~*'^'"^/L and summing over m and using (IB.27D we obtain the desired 
representation 



L 



n|j(|/)$(A) =« j dx{y,c\{ib[A'T^])-^\x,a)p\x)^{A) , (B.35) 


where 

{y,c\{ib[A^T,]r'\x,a) = T.'<MK].<A^) , = ^F^'^''^'^ • (^-36) 
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With Eq. (1B.35P one finds for the Coulomb Hamiltonian He defined by Eq. (I5.9P in this 
gauge with Jpp = Jd the following expression 

Hc = yI dxdyp\x)F^\x, y)p\y) (B.37) 
with the Coulomb kernel given by 

F''\x,y) = {x,a\{zb[A'T,]))-'\y,b) = ^'(x|n)(a|a)A;J(or|6)(n|y) . (B.38) 



Contrary to the Coulomb kernel in the pure Coulomb gauge (|B.19l) here only the zero 
mode m = 0" = is excluded, as indicated by the prime, while the mode m = 0, a = ±1 is 
included. The above considerations show that the Coulomb Hamiltonian depends on the 
details of the gauge fixing and is thus a priori not a physical quantity. 
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